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the  ocean,  received,  and  analayzed.  The  achieved  bit  error  rate  was  almost  the  same  as 
predicted  theoretically. 


I  INTRODUCTION 


I.l  General 

During  the  last  two  decades,  impressive  success  has  been  achieved  in  com¬ 
municating  and  transmitting  data  between  any  two  places  on  earth.  Great  success 
has  been  achieved  in  transmitting  and  receiving  high  quality  video  pictures  from  the 
moon  and  even  from  farther  in  the  universe,  but  much  less  attenuation  has  been  paid 
to  underwater  communication. 

The  major  reason  for  the  difficulties  and  the  poor  performance  in  underwater 
communication  is  due  to  the  fact  that  the  acoustic  channel  is  very  complex.  High 
attenuation,  strong  multipath,  and  fading  effects  cause  the  underwater  acoustic  chan¬ 
nel  to  be  almost  useless  as  a  channel  for  transmitting  acoustic  data. 

As  part  of  a  program  undertaken  at  the  Marine  Physical  Laboratory,  an 
underwater  communication  link  between  two  freely  driftingsensors  was  required. 
The  basic  requirements  of  the  link  were  the  following: 

a)  Operating  frequency  10-20  kHz. 

b)  Operating  range  2  -  5  km 

c)  Bit  rate  of  1  kbit/sec 

d)  Source  level  170dB/lpPa 

e)  Water  depth  4000  m 

d)  Receiver  depth  300  m 

e)  Transmitter  depth  500  m 

The  objective  of  this  research  was  to  develop  an  underwater  acoustic  link 
that  gives  reliable  performance  under  the  channel  constraints. 
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Chapter  1  of  this  thesis  gives  a  brief  review  of  past  and  present  efforts  in 
developing  underwater  communication  systems.  In  Chapter  2  we  summarize  the 
results  of  three  at-sea  experiments  which  were  conducted  during  1989-90.  By  using 
the  experimental  results,  a  simplified  model  of  the  underwater  acoustic  communica¬ 
tion  channel  is  obtained.  Chapter  3  deals  with  modulation  techniques  that  are 
efficient  in  a  strong  fading  environment.  It  is  shown  that  under  the  assumption  of 
slowly  time-varying,  frequency-selective  fading,  the  chirp  DPSK  modulation  tech¬ 
nique  is  a  superior  modulation  scheme  for  combating  the  fading  phenomena. 
Chapter  4  and  5  propose  bit  synchronization  and  equalization  methods  that  are  suit¬ 
able  to  the  assumed  channel.  Chapter  6  summarizes  the  results  of  at-sea  experiments 
in  which  the  proposed  chirp  DPSK  modulation  scheme  was  implemented. 

1.2  Underwater  acoustic  telemetry  -  overview 

Two  years  before  Columbus  discovered  America,  Leonardo  da  Vinci  wrote: 
If  you  cause  your  ship  to  stop,  and  place  the  head  of  a  long  tube  in  the  water  and 
place  the  outer  extremity  to  your  ear  you  will  hear  ships  at  a  great  distance  from 
you.  Until  at  the  turn  of  this  century,  essentially  no  improvement  was  made  in  the 
science  of  underwater  sound.  Scientific  research  using  diving  bells  and  the  beginning 
of  the  submarine  generation  motivated  researches  to  seek  better  ways  to  communi¬ 
cate  underwater. 

Although  underwater  echo  detection  schemes  emerged  prior  to  World  War  I, 
World  War  II  ushered  in  the  modem  age  of  underwater  sound  exploitation.  But  even 
then  the  water  medium  was  used  primarily  for  sound  navigation  and  ranging 
(SONAR). 
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In  the  past,  most  of  the  incentive  for  underwater  communication  has  come 
from  the  requirements  associated  with  submarines.  Now  there  are  very  rapidly  grow¬ 
ing  requirements  for  commercial  applications,  wireless  command  and  control  of 
unmaned  vehicles,  high  data  rate  video  transmission,  and  communication  with  deep 
divers  and  scientific  instruments  moored  deep  in  the  ocean. 

Conventional  telemetry  employs  either  electrical  signals  transmitted  over 
closed  wire  systems  or  some  form  of  electromagnetic  radiation  (radio  wave, 
microwaves,  etc.).  The  former  method  is  inconvenient  at  sea  and  latter  is  impractical 
due  to  the  high  attenuation  imposed  upon  the  entire  electromagnetic  spectrum  by  sea 
water.  Figure  1.1  5  shows  that  only  the  visible  region  holds  any  promise  for  underwa¬ 
ter  telemetry  and  even  then  in  the  practical  environments,  the  measured  values  of  rhe 
absorption  coefficients  far  exceed  the  clean  water  values. 
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Figure  1.1  Attenuation  of  electromagnetic  radiation  in  clear  water5. 
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Therefore,  it  is  clear  that  the  carrier  signal  must  be  an  acoustic  transmission 
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since  the  attenuation  imposed  by  sea  water  is  much  less  severe  provided  that  the  fre¬ 
quency  is  restricted  to  below  100  kHz  (See  Figure  1.2  [  Urick85  Fig  5.3  J) 


Figure  1.2  Low  frequency  acoustic  attenuation  coefficient85 

Although  prone  to  echo  interference,  signal  spreading  and  Doppler  frequency  shift, 
acoustic  signaling  is  the  oldest  and  the  most  successful  form  of  cableless  underwater 
communication,  because  of  the  relatively  low  attenuation  in  the  water. 

The  first  underwater  communication  system  was  a  telephone  developed  in 
1945  at  the  Naval  Underwater  Sound  Laboratory.  This  underwater  telephone  had 
been  designed  to  communicate  with  submerged  submarines  and  employed  the  upper 
sideband  of  an  8.3  kHz  U.S.S.B  system.  In  the  subsequent  15  years,  similar  systems 
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appeared  using  the  frequency  range  of  8-11  kHz  with  analog  modulation,  most  of 
them  using  AM  or  SSB  modulation.  Considerable  amount  of  transmitter  power  up  to 
some  hundreds  of  watts  was  employed  to  obtain  a  range  of  a  few  kilometers.  From 
these  first  operational  systems,  military  underwater  acoustic  communication  has 
moved  toward  lower  frequencies  that  permit  transmission  over  longer  ranges  or 
higher  frequencies  which  enable  higher  information  rate. 

Over  the  years,  various  forms  of  acoustic  communication  systems  have  been 
developed.  These  have  included  direct  AM  and  SSB  for  underwater  telephones,  or 
FM  for  sensor  data.  The  difficulties  of  AM  in  the  presence  of  multipath  and  disper¬ 
sion  were  recognized  and  experiments  with  frequency  modulation,  analog  as  well  as 
digital,  have  been  made.  Before  shifting  to  digital  modulation  techniques,  research¬ 
ers  have  partially  succeeded  in  transmitting  a  wide  variety  of  information.  A  direct- 
path  FM  acoustic  television  link  was  demonstrated  by  the  Ball  Brothers  Company  in 
1968,  where  an  information  rate  of  tens  of  kHz  was  achieved  for  several  kilometers 
of  water  path.  Also,  of  the  more  secure  command  control  links,  with  data  rates  as  low 
as  1  bit/sec  or  less  but  with  high  reliability  have  been  achieved  in  severe  multipath 
environment. 

At  the  same  time,  experiments  in  telemetry  over  a  vertical  path  from  sensors 
lowered  down  to  a  hydrophone  suspended  from  a  surface  ship  were  in  progress.  As 
the  phenomenon  of  multipath  is  less  dominant  in  the  vertical  direction,  reasonable 
results  have  been  achieved.  A  21  kHz  carrier  amplitude  modulation  was  used  in  a 
system  which  was  used  to  communicate  between  a  surface  ship  and  a  1800  foot  deep 
hydrophone29.  Tested  in  1966,  Performance  of  another  system,  using  FSK  modula¬ 
tion  with  a  carrier  frequency  of  40  kHz  and  achieving  data  rate  up  to  45  bit/sec  was 
reported  by  Hearn40. 
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Over  the  last  several  years,  there  have  been  many  attempts  at  designing  more 
reliable  communication  systems  using  digital  modulation  techniques.  It  became 
clear  that  analog  as  well  as  digital  modulation  techniques  using  amplitude  modula¬ 
tion  perform  inadequately  in  most  underwater  applications  because  of  the  destructive 
nature  of  the  underwater  channel.  Therefore  most  of  the  effort  has  been  directed 
toward  digital  systems  which  use  different  kinds  of  frequency  modulation. 

When  low  data  rate  was  required  for  control  telemetry,  undersea  navigation, 
and  positional  systems,  relatively  simple  underwater  acoustic  communication  sys¬ 
tems  yielded  reliable  performance.  Problems  arise  when  data  rates  increase  or 
transmission  have  long  range  multipath  echos.  The  destructive  nature  of  the  acoustic 
channel  causes  significant  interference  which  severly  limits  the  achievable  system 
performance.  In  the  past,  many  techniques  for  combating  the  multipath  problem 
have  been  proposed.  At  low  data  rates,  careful  selection  of  antimultipath  scheme  will 
suffice.  High  data  rate  telemetry,  however,  is  a  special  case  in  that  many  of  the  con¬ 
ventional  methods  cannot  be  employed. 

High  data  rates,  on  the  order  of  several  kbits/sec,  require  correspondingly 
high  transmission  bandwidths.  This  forces  the  use  of  high  carrier  frequency  signals. 
Unfortunately,  rapid  attenuation  of  the  acoustic  signal  limits  these  frequencies  to  no 
more  than  a  couple  of  hundred  of  kHz.  Over  the  last  15  years,  digital  techniques  for 
underwater  communication  have  been  widely  used  and  few  underwater  high  data  rate 
communication  systems  have  been  designed  and  developed  for  command  and  con¬ 
trol,  speech  transmission,  video  data  and  telemetry  information.  Pioneering  work  has 
been  done  at  the  Woods  Hole,  Oceonographic  Institute2, 7^ 

The  nature  of  the  underwater  acoustic  channel  prevents  the  use  of  coherent 
modulation.  Hence  non-coherent  FSK  modulation  is  the  most  commonly  used  modu¬ 
lation  method47.  The  MT-300,  designed  and  built  at  HonewelP6.  for  a  wireless  high 
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security  and  high  reliable  control  system,  achieved  good  performance  using  FSK 
modulation  with  frequency  diversity.  Tested  in  1982,  the  ATDL  79  is  another  acous¬ 
tic  telemetry  system  which  uses  noncoherent  FSK  modulation.  This  achieves  high 
reliability  by  using  a  smart  answer  back  communication  protocol. 

Significant  progress  in  underwater  communication  systems  is  achieved  when 
anti  multipath  modulation  and  coding  techniques  are  used.  In  1981,  Wax88  demon¬ 
strated  an  MFSK  system  using  a  15  kHz  carrier  with  a  bandwidth  of  3  kHz.  A  low 
data  rate  was  achieved  (40  bits/sec)  over  a  long  distance  (up  to  three  miles)  with  very 
low  bit  error  rate  (10“5).  A  method  proposed  by  Pieper  et  al71  is  now  under  develop¬ 
ment  at  Woods  Hole  15,16  where  a  data  rate  of  10  kbits/sec  over  10  km  in  shallow 
water  is  expected  to  be  achieved. 

Another  method  of  combating  the  multipath  problem  is  by  using  spread 
spectrum  techniques.  A  commercial  telemetry  system  using  frequency  hopping  has 
been  proposed  by  Datasonics.  A  different  spread  spectrum  technique  implementation 
for  underwater  communication  has  been  presented  by  W.  Hill.  42  This  system  uses  a 
chirp  FSK  modulation  as  a  way  to  combat  multipath  and  fading. 

Less  successful  results  have  been  achieved  in  developing  systems  for  acoustic 
telemetry  of  video  information.  Since  relatively  high  bit  rate  is  required  to  transmit 
video  information,  reasonable  performance  is  achieved  only  in  limited  scenarios. 

The  first  attempt  at  sending  TV  pictures  through  the  water  appears  to  be  the 
Cutlink  experiments10  performed  in  1968-9.  Both  analog  and  digital  FM  signals 
modulated  a  carrier  frequency  of  15  kHz  with  a  total  bandwidth  of  3  kHz.  Using  15 
watts  radiated  power,  a  data  rate  of  2.5  kbits/sec  for  a  vertical  path  of  up  to  6  km  was 
acieved.  Collins23  describes  an  acoustic  telemetry  system  for  video  information 
which  uses  a  narrow  beam  acoustic  link  to  transmit  slow  scan  video  frames  for  a  dis- 
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tance  of  up  to  30  meters.  The  transducer  projected  the  beam  at  an  angle  of  10°  down 
from,  horizontal  with  a  source  level  of  234.5  dB  @  1  p  Pa  at  lm.  In  the  presence  of 
weak  multipath,  good  quality  frames  were  received  continuously  at  8.5  seconds 
intervals. 

R.  M.  Dunbar30  reports  experiments  transmitting  T.  v'  pictures  at  very  high  data  rate. 
A  self-contained  TV  telemetry  system  which  works  in  the  vertical  direction  (from 
bottom  to  surface  and  vice  versa)  was  presented  by  B.  Leduc  and  G.  Ayela51  in 
1990.  This  system  works  at  21  kHz  and  uses  various  kinds  of  PSK  modulation.  The 
operating  range  is  up  to  6  km  depth  and  the  maximum  data  rate  is  9.6  kbit/sec  (19.2 
kbit/sec  at  2  km  depth).  The  effective  source  level  is  165  dB  @  1  Pa  at  1  m,  1  volt. 

Underwater  communication  equipment  has  a  wide  variety  of  applications 
and  many  companies  have  developed  specific  products  such  as  SONALINK50  and 
PASTY33. 


II  Channel  characterization  and  modeling 


II.I  General 


The  first  step  in  designing  a  complete  communication  system  is  to  character 
ize  the  medium  through  which  the  data  is  to  be  transmitted.  From  this  knowledge,  a 
model  of  a  communication  channel  can  be  developed  and  thus  an  expression  for  the 
received  signal  can  be  found 

The  ocean  is  far  from  being  the  ideal  medium  for  sound  propagation.  Its 
vertical  velocity  gradient  causes  the  transmitted  energy  to  be  refracted.  Some  times, 
communication  between  two,  points  can  be  impossible  due  to  the  formation  of  sha¬ 
dow  zones.  Also,  the  carrier  frequency  is  subject  to  a  randomly  fluctuating  Doppler 
shift  as  a  result  of  surface  and  internal  waves.  However,  the  signrl  is,  most  seriously 
distorted  by  multipath  propagation.  It  describes  a  condition  which  causes  the  signals 
from  a  given  transmitter  to  be  received  via  a  number  of  different  paths  with  each 
path  having  a  different  time  delay. 

Two  major  factors  contribute  to  this  multipath  propagation.  First  the  ocean 
is  inhomogeneous.  Its  own  wave  motion  breaks  up  the  surface  layer  and  produces  a 
turbulent  and  thermal  microstructure  consisting  of  a  random  distribution  of  small 
patches  of  water  with  slightly  differing  temperatures  which  are  continually  in 
motion.  This  inhomogenities  yield  variations  in  the  acoustic  refractive  index.  Thus  it 
is  possible  for  several  signals  each  radiated  in  slightly  different  initial  directions  to 
be  refracted  in  such  a  way  that  they  all  arrive  at  the  receiver  where  interference 
occurs  due  to  the  multiplicity  of  signal  transmit  times  and  hence  relative  phases 
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involved.  Often  refered  to  as  forward  scattering  or  fading.  This  phenomenon  pro¬ 
duces  random  fluctuations  in  both  the  amplitude  and  phase  of  the  carrier  signal 
which  become  progressively  more  troublesome  as  range  increases.  Usually,  the  sig¬ 
nal  amplitude  follows  a  Rayleigh  distribution.  In  some  special  cases  where  a  single 
multipath  is  particularly  strong,  the  signal  amplitude  follows  a  Rayleigh  Rice  distri¬ 
bution. 

The  second  factor  is  specular  reflection  from  boundaries  such  as  the  ocean 
surface,  bottom,  and  the  interface  between  the  water  and  man-made  objects.  Specular 
reflection  can  be  strong  approaching  the  strength  of  the  direct  signal  or  even  more55. 
The  transit  time  delays  associated  with  these  echoes  also  are  often  much  longer  than 
the  direct  signal.  Particularly  troublesome  are  the  reflections  from  the  surface  which 
not  only  fluctuate  greatly  but,  due  to  the  instantaneous  surface  shapes,  can  have 
amplitude  in  excess  of  the  direct  signal.  Multiple  reflections  also  are  possible  and 
usually  occur  between  the  surface  and  the  bottom  in  shallow  water.  Propagation 
delay  also  is  produced  by  reflection  from  the  boundaries  of  the  ocean.  Apart  from  the 
direct  path  it  is  also  possible  for  sound  to  travel  between  the  transmitter  and  the 
receiver  after  being  reflected  from  a  large  number  of  favorable  surfaces.  The  mul¬ 
tipath  phenomenon  results  in  the  reception  of  several  delayed  replicas  of  the  original 
transmitted  signal.  These  delayed  replicas  which  often  are  stronger  then  the  direct 
signal  itself  cause  strong  intersymbol  interference  (ISI). 

The  effect  of  multipath  propagation  on  the  signal  pulse  is  to  spread  it  in 
time,  the  pulse  spreading  being  determined  by  the  maximum  difference  in  the  transit 
delays  associated  with  the  various  paths.  Relative  movement  in  parts  of  the  medium 
which  causes  the  multipath  propagation  to  come  about  results  in  a  difference  in  the 
frequencies  of  the  signal  arriving  via  different  paths  (Doppier-shift).  Thus  a  modu¬ 
lated  signal  suffers  a  time  spreading  due  to  the  differential  time  delays  and  a 
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frequency  spreading  due  to  different  Doppler  shifts  arising  on  the  various  paths. 

As  mentioned  previously,  if  a  single  sine  wave  is  transmitted  through  such  a 
channel,  the  received  signal  is  the  resultant  of  a  number  of  sine  waves  of  slightly  dif¬ 
ferent  frequencies  which  have  been  subjected  to  different  time  delays.  As  a  result, 
the  amplitude  and  the  phase  of  the  received  signal  fluctuate  in  a  random  manner.  The 
amplitude  fluctuations  (fading)  of  the  sine  waves  of  different  frequencies  transmitted 
through  such  a  channel  are  correlated  if  the  frequency  difference  is  not  greater  than 
the  inverse  of  the  pulse  spreading  in  the  channel.  Hence,  a  band  limited  signal  with  a 
bandwidth  which  does  not  exceed  the  inverse  of  the  pulse-spreading  in  a  multipath 
channel,  retains  its  frequency  spectrum  (apart  from  an  ambiguity  arising  from  the 
Doppler  shift  and  the  Doppler  spreading),  when  transmitted  over  such  a  channel. 
Therefore,  the  inverse  of  the  pulse  spreading  can  be  termed  as  the  coherent 
bandwidth  of  the  channel.  If  the  signal  bandwidth  increases  beyond  this  value,  the 
various  frequency  components  of  the  signal  fade  in  an  uncorrelated  manner,  causing 
the  signal  to  be  distorted. 

If  all  the  spectral  components  of  the  signal  fluctuate  in  a  correlated  manner 
flat  fading  of  the  signal  occurs  and  the  channel  is  defined  as  a  frequency  nonselective 
channel.  Lack  of  correlation  of  the  fluctuation  characteristics  of  different  frequency 
components  signifies  frequency  selective  fading  of  the  signal.  The  underwater 
acoustic  channel  is  a  time  varying  channel.  Fading  and  multipath  behavior  change 
with  time.  A  channel  in  which  the  changes  are  very  slow  is  defined  as  a  slowly  time- 
varying  channel.  A  channel  which  changes  its  behavior  rapidly  is  defined  as  a  fast 
time  varying  channel. 

The  behavior  of  the  channel  in  the  time  domain  depends  on  the  operating 
scenario  and  the  environment  of  the  communication  system.  If  the  receiver  and  the 
transmitter  do  not  change  their  position  during  the  operation,  the  channel  is  a  very 
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slowly  time-varying  channel,  and  the  time  constant  of  the  channel  is  in  the  order  of 
some  hundreds  of  milliseconds.  If  the  receiver  and  /or  the  transmitter  move  and 
change  their  position  during  the  operation,  the  channel  becomes  a  fast  time-varying 
channel. 


II.2  Channel  characterization 


The  effect  of  fading  on  the  transmitted  signal  is  usually  modeled  as  a  linear 
time-varying  filter.  By  employing  a  narrowband  signal  model,  the  input  to  the  chan¬ 
nel  is  given  by 

S(O  =  Rejs(Oexp(0M)j  ■  (2.1) 

when  Re(  ■ )  denotes  the  real  part  of  ( • ) .  The  output  is  given  by: 

r(/)  =  Rej  f(Oexp(coeoj  (2.2) 

were: 


Y(t)  =  aS(t)+  J  h{t-v,t)S{x)dx  +  K(t)  ,  (2.3) 

and,  h{v,t)  is  the  lowpass  (complex)  impulse  response  of  the  fading  channel,  x 
represents  the  usual  filter  response  variable  and  the  t  -  dependence  indicates  that  the 
very  structure  of  the  impulse  response  changes  with  time84.  h{x\t)  is  modeled  as  a 
complex,  zero  mean  Gaussian  random  process.  Hence  it  has  Rayleigh  amplitude  and 
uniform  phase  distribution.  The  Fourier  transform  of  h(x;t)  with  respect  tot 
is  the  time- varying  equivalent  lowpass  transfer  function  of  the  channel.  H(f\t)  is  also 
a  zero  mean,  complex  Gaussian  random  process.  In  addition  to  the  transfer  function, 
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the  channel  (or  the  fading  process)  is  characterized  by  the  ensemble  autocovariance, 
RT{zui2,t\ti),  of  the  channel  impulse  response,  h(x,t).  In  most  cases,  the  channel  is 
assumed  to  be  wide  sense  stationary  in  both  dimensions,  i.e. 

R~(x,  =  RT(z:At) 

where,  t  =1^2  and  a t=tx-t2. 

Another  expression  which  is  used  to  characterize  the  fading  channel  is  the 
spaced-tone  complex  covariance  which  is  defined  as: 

=  •  (2-4) 

This  expression  characterizes  the  relationships  of  values  of  the  transfer  function  at 
different  frequencies  and  at  different  times.  Under  the  assumption  of  stationary, 

A  measure  of  the  degree  of  frequency  selectivity  commonly  employed  in  the  litera¬ 
ture  is  the  distance  pH  between  the  j  points  of  the  frequency  correlation  function35. 

In  most  cases, 

R/r (A/  ;At )  =  Rf  (A/ )  Rt(At )  .  (2.5) 

For  R,(At)  which  is  almost  constant  during  the  observation  time,  the  fading  is  known 

as  very  slowly  time  varying  fading.  When  R/(Af)  is  constant  (i.e  not  frequency 
dependent),  the  fading  is  known  as  frequency  nonselective  fading.  When  Rf  (A/ )  is 
such  that  pH  is  much  smaller  than  the  data  bandwidth,  the  fading  is  known  as  fre¬ 
quency  selective  fading.  The  ratio  between  the  available  channel  bandwidth  and  pH 
is  the  amount  of  diversity  that  the  channel  provides  and  is  defined  as  K .  In  the  most 
general  case  of  fading,  the  received  signal  Y(t )  consists  of  two  components:  a  single 
specular  component  and  a  diffuse  Rayleigh  distributed  component.  When  a  in  (2.3) 
is  zero,  only  the  diffuse  component  of  fading  exists.  In  this  case,  Y  in  (2.4)  is  a  Ray¬ 
leigh  random  process  with  the  following  distribution  function: 
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P(y)  =  -^;  exp 


2  a2 


U(y), 


(2.6) 


where  U(  ■ )  is  the  unit  step  function.  Defining  z  =  y  ,  we  have 


P(j)=  A-  exp(-z)  U(z)  .  (2.7) 

The  fading  phenomena  has  been  well  explored  in  connection  with  elec¬ 
tromagnetic  communication  links.  Most  of  the  effort  has  been  directed  toward  a  deep 
understanding  of  the  fading  phenomena  and  its  aspects  in  frequency  ranges  used  in 
radio  communication.  Relatively  little  information  has  been  published  in  the  litera¬ 
ture  describing  and  analyzing  the  fading  phenomena  and  the  character  of  the  under¬ 
water  acoustic  communication  channel. 


II.3  Channel  character  -  experimental  results 


During  1989-1990,  three  experiments  were  conducted  at  sea.  The  objectives 
of  the  experiments  were  *o  measure  the  transmission  characteristics  of  the  acoustic 
channel  at  high  frequency  (10  to  20  kHz)  and  to  explore  the  nature  of  fading  and 
multipath  in  shallow  water,  midrange  acoustic  channel.  The  first  two  experiments 
were  carried  out  on  March  16,  1989  and  on  December  20,  1989  and  were  located  at 
32°40'N  117°35.6'  W.  The  third  experiment  was  conducted  on  May  22,  1990  and  was 
located  at  32°50'N  117°35'W. 

During  the  three  experiments  the  sea  state  was  between  zero  and  one  and  the 
wind  speed  was  between  5  and  12  knots.  The  purpose  of  each  experiment  was  to 
transmit  a  set  of  waveforms  from  a  transducer  deployed  deep  in  the  ocean  from  a 
ship  and  receive 
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Figure  2.1  Experiment  setup 


Figure  2.2  Projector  and  Monitor  Hydrophone  Location 
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the  transmitted  signal  with  four  sonobuoys  located  1  km  apart  and  transmitting  the 
received  signal  via  an  RF  link  back  to  the  ship  (see  Figure  2.1  and  2.2).  The 
transmitted  and  the  received  signals  (from  the  sonobuoys)  were  recorded  simul- 
tanously.  In  addition,  a  monitor  hydrophone  was  deployed  close  to  the  projector  and 
provided  a  replica  of  the  waveform  which  was  transmitted  through  the  water. 

Among  the  many  waveforms  transmitted,  two  specific  waveforms  enabled  us  to 
measure  the  channel  charactristics: 

(1)  A  set  of  nine  tones  equally  spaced  was  transmitted  for  2  minutes.  This  set 
of  tones  was  retransmitted  five  times,  where  each  time  the  space  between  the 
tones  was  changed  (50,  100,  250,  500  and  1000  Hz  appart).  The  fifth  tone 
was  always  centered  at  15  kHz. 

(2)  A  set  of  spread  spectrum  waveforms  implemented  by  chirp  signals  was 
transmitted.  Five  different  waveforms  were  transmitted,  1  msec  chirp  with 
time-bandwidth  product  of  10,  2  msec  chirp  with  time  bandwidth  of  20,  4 
msec  chirp  with  time  bandwidth  of  40,  8  msec  chirp  with  time  bandwidth  of 
80,  and  16  msec  chirp  with  time  bandwidth  of  160.  All  these  waveforms 
were  designed  such  that  they  occupied  a  bandwidth  of  10  kHz  between  10 
and  20  kHz.  Figure  2.3  shows  the  8  msec  chirp  and  its  power  spectrum. 
Each  waveform  was  transmitted  every  four  seconds  for  2  minutes. 
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Figure  2.3  Time  series  and  power  spectrum  of  an  8  msec  long  chirp 
waveform. 


The  channel  parameters  measured  were:  (1)  the  channel  multipath  character,  (2)  the 
coherence  function  as  a  function  of  time,  (indicates  the  nature  of  the  time  varying 
character  of  the  acoustic  channel),  and  (3)  histograms  of  amplitude  variation  at  each 
frequency  commponetu  of  the  waveforms  (enable  us  to  measure  the  statistics  of  the 
fading  and  the  correlation  between  two  frequency  componnents). 
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II.3.1  The  channel  multipath  characteristcs 

The  waveforms  that  are  received  by  the  sonobouys  are  distorted  replicas  of 
the  transmitted  signal  plus  many  other  delayed  replicas  that  result  from  the  mul¬ 
tipaths  (see  Figure  2.4).  The  multipath  characteristics  of  the  channel  are  calculated 
in  two  different  ways.  The  first  is  based  on  inverse  Fourier-transforming  the  ratio  of 
the  received  and  the  transmitted  chirp  waveforms  spectra.  The  second  is  by  correlat¬ 
ing  the  received  and  the  transmitted  waveforms  and  envelope-detecting  the  result. 
Typically,  the  chanel  multipath  characteristcs  are  caracterized  by  the  received  signal 
from  the  direct  path  being  much  weaker  than  the  signal  received  from  the  first  mul¬ 
tipath  (caused  by  reflection  from  the  sea  surface)  (see  Figure  2.5).  The  other  mul¬ 
tipaths  having  larger  delay  are  caused  by  reflection  from  the  sea  bottom  or  from  mul¬ 
tiple  reflections  (bottom  and  surface)  and  are  much  weaker  than  the  first  multipath. 


Normalized  magnitude  (dB) 
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Figure  2.4  Typical  received  signal 
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Figure  2.5  Typical  channel  multipath  character 
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The  underwater  acoustic  channel  is  a  time-varying  channel  and  hence,  the 
channel  multipath  character  also  varies  with  time.  Figures  2.6  and  2.7  present  the 
multipath  character  of  the  channel  as  calculated  from  twenty  eight  consecutive  eight 
msec  chirp  signals  transmitted  every  4  seconds.  The  first  20  msec  (the  main  path  and 
two  first  multipaths)  of  each  multipath  character  is  shown  in  Figure  2.6.  The  channel 
multipath  character  in  the  time  interval  between  370  msec  and  530  msec  is  shown  in 
Figure  2.7.  These  two  plots  (which  are  typical  examples  of  the  three  at-sea  experi¬ 
ments)  illustrate  the  time  variation  of  the  channel  multipath  character.  Careful 
inspection  of  the  first  20  msec  (Figure  2.6)  shows  not  only  that  the  multipath  inten¬ 
sity  changes  with  time  but  that  the  multipaths  die  and  rebuild  during  a  long  time  of 
inspection.  Another  interesting  property  of  the  channel  multipath  character  is  that 
each  multipath  consits  of  a  group  of  many  close  multipaths  (see  Figures  2.5  and  2.6) 
(micromultipath  effect).  This  effect  is  more  dominant,  as  the  number  of  reflections 
(from  bottom  as  well  as  from  surface)  become  larger  (see  Figure  2.7). 
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Figure  2.7  Typical  time  variaton  of  channel  multipath  character, 
(time  slice  between  370  and  530  msec) 
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II.3.2  Time  variation  properties  of  the  channel 

One  of  the  characteristics  of  an  underwater  acoustic  channel  is  its  time  vary¬ 
ing  property.  The  squared  magnitude  coherence  (SMC)  is  a  good  measure  of  the 
correlation  between  the  channel  character  in  two  different  time  intervals28.  A  spread 
spectrum  waveform  between  10  and  20  kHz  (8  msec  chirp)  was  transmitted  every  4 
seconds  and  the  squared  magnitude  coherence  function  between  the  first  received 
waveform  and  the  rest  of  the  waveforms  was  calculated  (Figure  2.8).  From  the  calcu¬ 
lated  SMC  it  is  observed  that  for  most  frequency  components,  the  channel  properties 
are  rearly  time  invariant  and  the  coherence  is  very  high  (almost  1).  An  exception  is 
in  the  frequency  range  between  13  and  16  kHz  (normalized  frequencies  of  0.26  and 
0.32),  where  the  coherence  becomes  smaller  as  the  time  difference  becomes  larger. 
As  a  result,  one  can  conclude  that  the  rate  of  variation  of  the  channel  character  is 
much  lower  than  the  bit  rate  (seconds  compared  with  miliseconds)  and  that  the  chan¬ 
nel  can  be  assumed  to  be  very  slowly  time-varyingl.  A  different  way  to  show  the 
rate  of  variation  of  the  channel  character  as  the  function  of  time  is  to  look  at  the 
spectrum  variation  of  the  received  signal  as  a  function  of  time  (see  Figure  2.9).  The 
spectrum  has  small  variation  in  a  time  interval  of  some  tens  of  seconds  which  showse 
again  that  the  channel  is  very  slow  time-varying. 
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Figure  2.8  The  square  magnitude  coherence  function  of  the  channel 
as  a  function  of  time. 
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II.3.3  The  frequency  selectivity  properties  of  the  channel 

Another  characteristic  featuic  of  a  communication  channel  is  its  frequency  selec¬ 
tivity.  This  parameter  was  measured  by  transmitting  a  set  of  tones  and  calculating  the 
correlation  between  the  envelopes  of  the  tones. 

The  envelope  variation  of  each  frequency  component  is  presented  in  Figures  2.10 
and  2.11.  Figure  2.10  presents  the  envelope  variation  of  each  frequency  component 
for  tones  seperated  50  and  100  Hz  apart  and  Figure  2.1 1  for  250  and  500  Hz  appart. 
Figure  2.12  gives  the  envelope  variation  when  the  tones  are  seperated  1000  Hz  apart. 
The  correlation  coefficient  of  the  envelope  variation  between  the  nine  tone  groups 
was  calculated  with  respect  to  the  center  tone  (15  kHz)  (see  Figure  2.13).  The  sta¬ 
tistical  distribution  of  the  amplitude  of  each  tone  was  calculated.  Figures  2.14  and 

2.15  present  the  histograms  of  the  envelope  of  the  tones  when  the  tones  are  250  Hz 
and  1000  Hz  apart. 

Inspection  of  the  figures  shows  that  the  correlation  of  the  fading  phenomena 
between  two  frequencies  is  time  dependent.  The  correlation  between  the  same  two 
frequencies  obtains  different  values  at  different  time  intervals.  As  an  example,  the 
correlation  between  signals  at  15  and  15.1  kHz  is  -0.039  in  one  case  (tones  50  Hz 
apart)  and  0.1111  in  the  second  (tones  100  Hz  apart).  The  correlation  between  two 
tones  500  Hz  apart  (15  and  15.5  kHz)  is  0.15564  in  one  time  interval  and  -0.0452  in 
other  time  interval.  A  second  feature  of  the  data  is  that  not  only  the  correlation  func¬ 
tion  is  time-varying,  but  the  statistics  of  the  envelope  are  also  time-varying.  Figure 

2.16  gives  the  histogram  of  the  envelope  of  the  central  tone  (at  15  kHz)  at  five  dif¬ 
ferent  time  intervals  of  65  seconds  length  and  120  seconds  apart. 
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Figure  2.13  The  correlation  between  the  envelope  of  nine  tones  and  the  central  tone 
for  the  five  sets  of  tones. 
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Figure  2.14  Histograms  of  the  envelopes  of  the  250  Hz  spaced  tones. 
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Figure  2.15  Histograms  of  the  envelopes  of  the  1  kHz  spaced  tones 
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Figure  2.16  The  histograms  of  the  15  kHz  tone,  at  different  time  periods 
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From  Figure  2.16,  one  can  see  that  at  each  different  time  interval  the  statistics  of  the 
envelope  change  between  a  Rayleigh-like  distribution  and  a  Rayleigh-Rice  like  dis¬ 
tribution86.  The  Rayleigh-like  distribution  exists  when  no  specular  multipath  com¬ 
ponent  exists,  and  the  Rayleigh-Rice  like  distribution  appears  when  a  strong  specul- 
lar  multipath  component  exists. 

II.4  Channel  modeling 

The  analysis  of  the  three  at-sea  experiments  shows  that  the  acoustic  com¬ 
munication  channel  is  a  time-varying  frequency  selective  channel,  which  suffers 
from  strong  multipath  effects.  A  simplified  model  of  this  complicated  channel  is  a 
sum  of  independent  channels  each  having  a  different  transfer  function  and  a  different 
delay  depending  on  the  delay  caused  by  the  multipath  (see  Figure  2.17).  The  chan¬ 
nel  is  assumed  to  be  a  frequency  selective  Rayleigh  fading  channel  with  low  correla¬ 
tion  between  two  frequencies  if  they  differ  by  50-100  Hz. 


Figure  2.17  A  simplified  model  of  an  underwater  acoustic  communication  channel. 
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The  equivalent  lowpass  multipath  character  of  the  channel  can  be  described 
as 

hc  (t.0  =  hQ(u )  +  £  ^  (r -TD,  j )  (2.8) 

where 

f/.-( w)  t«  TD, 

/i0(t,r)  is  the  direct  path  and  h, (x,t ),  i  =  are  the  delayed  (with  delay  of  TD,-) 

pathss  caused  by  the  multipaths. 

A  more  simplified  model  for  the  channel  can  be  described  as 

hc  (x,f )  =  Oo  6(  x )  +  £  a,  z  -  TD, ).  (2.9) 

i=  I 

As  seen  from  Figures  2.4,  2.5  and  2.6,  one  of  the  multipath  signal  is  much  stronger 
than  the  direct  path.  This  result  was  also  verified  experimantally  and  explained  by 
Lord  and  Plemons  (see  their  Figures).  55  For  cases  where  an  a,  exsits  such  that 
a,  >ao,  the  transfer  function  of  the  channel  is  a  non  minimum  phase  transfer  function. 

m 

A  necessary  condition  for  a  function  to  be  a  minimum  phase  is  that  £  |  a,  |  <  |  oto  I  . 

1=1 
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III.l  Introduction 

The  underwater  acoustic  communication  channel  has  been  described  and 
modeled  in  Chapter  2.  As  fading  is  a  dominant  phenomena,  when  the  modulation 
approach  is  selected,  the  ability  to  combat  destructive  fading  must  be  considered. 
The  stochastic  nature  of  the  channel  destroys  any  phase  coherence  and  therefore  only 
incoherent  modulation  can  be  considered. 

The  simplest  way  to  combat  the  fading  channel  is,  to  select  signal  parame¬ 
ters  (symbol  bandwidth,  keying  rate)  that  allow  essentially  distortionle.  s  reception  of 
individual  symbol  waveforms.  This  is  possible  for  unspread  channels,  by  selecting  a 
symbol  waveform  with  roughly  unity  time-bandwidth  product  and  for  which  the 
symbol  duration  T  satisfies  84 

1  -.  1 
—  «  T  «  — — 

Ph  bd 

Pi,  is  defined  in  Chapter  II. 2  and  B0  is  the  Doppler  spread  of  the  channel.  For  a 
keyed  stream  that  does  not  satisfy  the  inequality  above,  other  methods  such  as  diver¬ 
sity,  adaptive  equalization,  coding,  and  special  waveforms  and  modulation  must  be 
used. 

III.2  Frequency  diversity  techniques  for  a  selective  fading  channel 

Almost  all  diversity  techniques  are  based  on  the  fact  that  errors  occur  when 
the  channel  is  in  a  deep  fade.  Supplying  the  receiver  with  several  replicas  of  the 
same  information  signal  transmitted  over  independently  fading  channels  will  reduce 
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the  probability  that  all  the  signals  will  fade  simultaneously. 

The  main  diversity  methods  are  75,81,84  :  (1)  space  (spaced  antenna)  diversity,  (2) 
angle  (of  arrival)  diversity,  (3)  polarization  diversity  (not  relevant  in  underwater 
acoustic  communication),  (4)  frequency  diversity  and  (5)  time  (signal-repetition  (bit 
interleaving))  diversity. 

A  detailed  analysis  of  several  frequency  diversity  methods  based  on  coherent 
and  differential  PSK  and  coherent  and  incoherent  FSK  modulation  is  given  by 
Proakis  75  (Chapter  7  )  and  by  Stein**1.  Another  frequency  diversity  method  used 
successfully  in  radio  communication  fading  channels  is  spread  spectrum  direct 
sequence  modulation.  However,  the  need  for  sequence  acquisition  and  maintenance 
of  synchronization  while  one  of  the  terminals  is  in  motion  causes  serious  problems 
for  this  technique  in  an  underwater  communication  system. 

In  this  section,  we  will  summarize  the  performance  of  the  most  common  fre¬ 
quency  diversity  methods  in  the  presence  of  frequency  selective,  very  slowly  time- 
varying  Rayleigh  fading. 

II  1.2.1  The  combining  method 

Figures  3.1  and  3.2  give  block  diagrams  of  coherent  and  incoherent  binary 
FSK  receivers  with  combined  L  branches.  The  receivers  are  built  from  L  duplicates 
of  the  same  branch  tuned  to  the  L  different  carrier  frequencies.  Each  branch  is 
matched  to  the  appropriate  transmitted  signal. 
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It  is  easy  to  show  that  when  each  branch  independently  fades,  the  optimal  maximum 
likelihood  decision  rule  is  to  sum  the  outputs  of  all  the  mark  channels  and  all  the 
space  channels  of  the  branch  receivers  and  to  compare  them.  Choosing  the  mark  and 
the  space  at  each  branch  such  that  they  will  fade  together  will  significantly  simplify 
the  optimal  receiver. 

For  a  BPSK  receiver  with  combined  L-th  order  diversity  and  a  very  slowly 
time-varying  frequency  selective  Rayleigh  fading  channel,  the  bit  error  probability 
is75  [  p.  723  ] 


P. 


(3.1) 


where 


Ye 


is  the  average  signal  to  noise  ratio  for  all  channels. 


When  a  very  slowly  selective  fading  is  assumed,  the  instanstenous  signal  to 
noise  ratio  at  each  branch  remains  constant  for  a  long  period  of  time  relative  to  the 
message  length.  In  such  a  case,  yc  is  considered  as  the  average  signal  to  noise  ratio 
over  all  branches  and  since  the  transmitted  power  is  divided  among  all  L  branches, 
we  have: 


i  =  lZ*  =  T  (3-2) 

where  y,  denote  the  signal  to  noise  ratio  at  the  /th  branch  of  the  receiver. 

Under  the  assumption  that  the  mark  and  the  space  frequencies  are  close 
enough  that  they  fade  together  and,  again  assuming  a  very  slowly  time-varying  fre¬ 
quency  selective  Rayleigh  fading,  we  find  that  the  bit  error  probability  is  the  same  as 
(3.1)  (BPSK)  with  q  as  follows: 


(a)  foi  the  coherent  case 
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Figure  3.4  Bit  error  probability  of  binary  incoherent  FSK  with  L  order 
diversity  as  a  function  of  SNR 
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Figure  3.5  Bit  error  probability  of  binary  coherent  FSK  with  L  order  diversity  as  a 
function  of  SNR. 


IIL2.2.  Switch  diversity 

Switch  diversity  is  a  method  where  only  the  branch  having  the  largest  signal 
to  noise  ratio  72  is  used  at  any  given  instant  for  the  decision  rule. 

For  a  Rayleigh  fading  channel. 


P(y,)=  J-expJ-24  (3.5) 

yc  l  YcJ 

If  we  define  y  as  the  maximum  from  all  y,,  then  for  a  switch  diversity  receiver  with  L 


branches. 


p  (y)  =  -r-  exp  (-4-)  | \-exp 


•3 


and 


Pt  =  jpt/J.e)py(y)dy. 


(3.7) 


For  binary  PSK  with  an  L  branch  switch  diversity  receiver,  the  bit  error  probability  is 


£-1 


P.  =  f  exp(-^-) \  l-exp(-^-)l  <s>(-J2y)dy 

0  yc  7c  I  7c  I 


(3.8) 


where 


<K*)  =  |  zxy(-Ay2)dy 


and 


-\2y 


Pt  =  J  £  ^  (“O*  exP  (--^)  J  exp  (’/iy2)  dy  rfy  •  (3-9) 

wyc  o  t=o  (  *  J  7c  -i 

Changing  the  order  of  summation  and  integration,  and  integrating  by  parts,  we  obtain 


(3.10) 


For  coherent  binary  FSK  modulation,  substituting  'Ay  instead  of  y  in  (3.8),  will  give 
us  the  following  bit  error 


£.-1 


P.  =  f  I  *  (-D* 

z  t=0 


7.-1) 

k 


-v 


27+yc 


For  incoherent  FSK 


(3.11) 


p.  =  I  exp(-^r-)'|  1-exp  (-^-) 
27c  o  7c  I  7c 


L-l 


exp(-y)  £*7 


(3.12) 


27c  *=o  (  *  J 


2+7c+2  kj 
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III.3  Special  modulation  techniques 


m.3.1  Combined  modulation  and  coding 

Figure  3.6  shows  a  block  diagram  that  describes  the  idea  of  combining  cod¬ 
ing  and  modulation71. 


Input  Data 
Sequence 

0011001001 


Coder 


oi  i. ./i  oo.  y 


Serial 

to 

Parallel 


0 

4- 

1 


F 


Modulator 


Figure  3.6  Combining  coding  and  modulation  -  block  diagram. 


A  stream  of  binary  data  bits  is  supplied  to  the  encoder  which  encodes  the 
information  bits  (most  often  to  a  (n,k)  block  code).  The  modulator  accepts  a  block  of 
n  bits  corresponding  to  a  codeword  and  assigns  each  bit  to  a  cell  in  the  partitioned 
signal  space.  Waveforms  are  constructed  by  following  the  convention  that  a  tone 
pulse  is  generated  in  a  cell  if  a  "one  "  is  assigned  to  that  cell  and  that  no  energy  is 
transmitted  in  a  cell  to  which  a  "zero"  is  assigned.  This  choice  of  basic  cellular 
modulation  is  based  on  our  assumed  inability  to  detect  received  tones  coherently. 
The  waveform  type  that  results  for  an  ensemble  of  such  signals  is  termed  multitone 
on/off  keying  (MTOOK).  An  alternative  mapping  of  the  n  bits  into  channel 
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waveforms  can  be  accomplished  by  means  of  frequency  shift  keying.  With  FSK, 
each  bit  in  a  codeword  is  assigned  two  cells:  one  cell  for  a  "one"  and  the  other  for  a 
"zero". 

The  coded  waveform  technique  has  a  major  advantage  over  the  regular  fre¬ 
quency  diversity  methods  because  it  combines  coding  gain  and  frequency  diversity 

^min 

together.  The  obtained  diversity  order  (L)  is  — —  where  d^m  is  the  minimum  Ham¬ 
ming  distance  of  the  code. 

III.3.2  Chirp  as  a  frequency  diversity  technique 

At  least  25  years  ago93,  it  already  was  recognized  that  a  modulated  chirp  sig¬ 
nal  would  be  a  simple  and  effective  tool  for  combating  the  problem  of  fading,  mul¬ 
tipath,  5,6,37  impulse  noise44,87,  and  cross-talk25.  Furthermore,  some  successful 
experiments  were  carried  out  implementing  chirp  signals  in  underwater  acoustic 
telemetry42.  It  is  well  known  5,84  that  a  communication  system  based  on  coherent 
signal  detection  is  useless  in  a  random  fading  environment.  Thus,  among  all  the 
modulation  techniques  based  on  chirps,  only  methods  based  on  incoherent  modula¬ 
tion  are  used.  In  a  binary  communication  system,  the  mark  and  space  are  defined  by 
two  different  signals.  In  a  chirp  modulation  system,  these  two  signals  can  occupy 
different  frequency  bands  with  the  same  or  inverted  slopes  (frequency  vs.  time  char¬ 
acter),  or  can  occupy  the  same  frequency  band  but  with  inverted  slopes.  Bemi  & 
Gregg6  found  the  optimal  signal  sets  for  the  case  where  the  two  signals  had  oppo¬ 
site  slopes  and  occupy  partially  overlapped  frequency  regions.  Zaytsev  & 
Zhuravlev96  analyzed  the  performance  of  a  FSK  chirp  modulation  system  where 
both  signals  (mark  and  space)  occupied  the  same  frequency  bandwidth  with  opposite 
slopes. 
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Since  a  very  slowly  time-varying  frequency  selective  Rayleigh  fading  chan¬ 
nel  has  been  assumed  and  the  relative  drift  between  the  transmitter  and  the  receiver 
is  assumed  to  be  very  small,  the  difference  between  received  phases  in  consecutive 
bits  when  coherent  modulation  is  used  is  negligible.  Therefore,  DPSK  chirp  modu¬ 
lation  /demodulation  can  be  used.  Figure  3.7  and  Figure  3.8  show  two  applications 
of  a  binary  DPSK  receiver.  Figure  3.7  is  the  most  common  and  the  most  often  used 
scheme. 


H  i 


Wo 


Figure  3.7  Suboptimum  binary  DPSK  demodulator. 


COS((Oq/  +  Vi\xi2) 


sin(aty  +  Vip/2) 

Figure  3.8  Optimun  binary  DPSK  demodulator. 


Although  it  has  been  shown  that  when  the  effect  of  the  bandpass  filter  is  taken  into 
consideration,  the  demodulator  shown  in  Figure  3.7  is  not  optimal.  Due  to  its  simpli¬ 
city  it  is  prefered  over  the  configuration  shown  in  Figure  3.8.  J.H.  Park69  proved 
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that  when  the  bandpass  filter  effect  is  considered,  the  demodulator  shown  in  Figure 
3.7  has  worse  performance  (about  1  to  2  dB)  than  the  optimum  receiver.  He  also 
showed  that  the  demodulator  shown  in  Figure  3.8  is  optimum.  When  the  effect  of  the 
bandpass  filter  is  neglected,  both  receivers  give  the  same  performance. 

In  Appendix  A,  we  show  that  in  the  presence  of  a  frequency  selective  Ray¬ 
leigh  fading  channel  with  reasonable  degree  of  diversity,  K ,  the  probability  of  error 
when  chirp  DPSK  modulation  is  used  turns  out  to  be  the  same  as  regular  DPSK 
modulation  in  the  presence  of  AWGN.  The  bit  error  probability  is  inversely  propor¬ 
tional  to  the  exponent  of  the  signal  to  noise  ratio  instead  of  being  proportional  to  the 
inverse  of  the  signal  to  noise  ratio  as  in  BPSK  in  the  presence  of  Rayleigh  fading. 
Binary  DPSK  chirp  modulation  as  well  as  higher  order  DPSK  can  be  used.  Figure  3.9 
shows  the  achieved  bit  error  probability  as  a  function  of  the  signal  to  noise  ratio54  for 
2,4,8,16  and  32-ary  DPSK. 
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Figure  3.9  Bit  error  rate  of  2,4,8,16  and  32-ary  DPSK54 

Let  us  define  the  bandwidth  expansion  of  a  modulation  techniqe  Be  as  the 
time-bandwidth  product  per  information  bit.  It  is  a  measure  of  the  efficiency  of  the 
modulation  technique.  As  explained  in  Appendix  A,  when  a  chirp  signal  is  used,  a 
time  bandwidth  product  greater  than  10  is  recommended  in  order  to  prevent 
significant  intersymbol  interference.  When  a  4-ary  chirp  DPSK  is  used,  the 
bandwidth  expansion  of  the  system  is  5  (for  BT=10).  Figure  3.10  gives  a  comparison 
between  different  types  of  waveforms  combined  with  coding  with  Be  =4  and  a  4-ary 
chirp  DPSK  ( Bt  =  5)  the  presense  of  Rayleigh  frequency  selective  channel. 
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Figure  3.10  Comparison  of  4-ary  chirp  DPSK  with  various  types  of  combined  modulated 

and  coded  waveforms75  (Chapter  7) 

Curve  1  shows  the  probability  of  bit  error  for  a  4-ary  chirp  DPSK  modulation 

scheme.  Curve  2  corresponds  to  a  rate  y  binary  convolutional  code  with  constraint 

length  5  and  with  soft  decision  decoding.  Curve  3  depicts  the  performance  of  a  rate 

y  dual  -2  with  block  orthogonal  [0(4,2)]  inner  code,  curve  4  corresponds  to  (20,5) 

combined  Hadamard  code  and  curve  5  coresponds  to  a  (24,12)  Golay  code  with  soft 
decision  decoding. 

From  Figure  3.10,  one  can  conclude  that  a  4  -ary  chirp  DPSK  modulation  is  superior 
to  the  combined  modulated  and  coded  waveform  with  similar  bandwidth  expansion. 
Moreover,  inspection  of  Figuie  A.  10  (Appendix  A)  shows  that  the  chirp  DPSK  also 
is  superior  to  the  BFSK  and  BPSK  modulation  methods  even  when  a  high  order  of 
frequency  diversity  is  used. 
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III.3.3  Unequal  bit  error  rate  method 

Under  given  constraints,  the  main  purpose  of  a  communication  system  is,  to 
transfer  information  between  two  points  with  as  high  a  reliability  as  possible.  In  con¬ 
ventional  communication,  analog  as  well  as  digital,  a  lot  of  effort  is  directed  toward 
methods  of  improving  the  reliability.  In  digital  communication,  sophisticated  modu¬ 
lation  schemes,  channel  equalization,  and  error  correction  are  used  to  improve  the 
system  performance. 

In  conventional  digital  communication  systems,  all  cf  the  information  sym¬ 
bols  of  a  message  are  regarded  equally  significant  and  equal  apriori  symbol  probabil¬ 
ity  is  assumed.  In  such  cases,  a  common  measure  of  the  system  performance  is  the 
achieved  symbol  (or  bit)  error  rate  at  the  receiver’s  output.  However,  in  some  cases, 
not  all  data  are  of  equal  value  for  the  user.  For  example,  an  analog  data  which  has 
been  sampled  by  an  A/D  convener.  Errors  in  the  higher  order  bits  (MSB)  of  numbers 
would  be  more  costly  than  errors  in  the  lower  order  bits.  In  such  cases,  when  each  bit 
has  a  different  significance,  the  symbol  (or  the  bit)  error  rate  is  not  a  good  measure  of 
the  system  performance  any  more  and  another  measure  has  to  be  defined. 

The  following  model  of  a  digital  communication  system  is  considered.  An 
analog  signal  is  sampled  and  transmitted  digitally  via  a  noisy  channel.  The  receiver, 
receives  the  digital  data  and  reconstructs  the  analog  data. 

The  sampled  signal  ( V(n ))  at  each  sampling  time  can  be  described  as 

^)4lc«2'  (3.13) 

*  i=0 

where  C„  is  either  one  or  zero. 

Each  sampled  signal  is  uniquely  described  by  its  b  bits  Cm  i  =0,1,2..../? -l , 
and  the  only  data  that  is  transmitted  are  the  C ’s.  Each  set  of  b  C,  ’s  that  describes  one 
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sample  of  the  sampled  signal  is  defined  as  a  word  (W).  The  receiver  estimates  the 
Cm ,  and  reconstructs  the  signal  by 

V(»)=-^-  I  C»2‘  ,  (3.14) 

where  Cm  is  the  estimated  value  of  C*  . 

Instead  of  using  the  bit  error  rate  as  a  performance  measure,  we  define  the 
mean  square  error  (m.s.e)  (e2)  as  a  new  measure  of  the  system  performance,  where 


e2H£|[V'(n)-y(n)]2j 
(£[•]  denote  expected  value  of  [•]). 

In  Appendix  C,  we  show  that 


^r£|[I(Cm-Cw)]22* 


(3.15) 


i=0 


e2=i  I.  Pi  2*  I3-16) 

1  i=0 

where  p,  is  the  probability  of  incorrect  detection  of  the  i-th  bit  in  the  word  (W).  In 
the  most  general  case,  the  bit  error  probability,  a,  at  each  bit  may  have  different 
values.  In  a  particular  case,  where  all  bits  have  the  same  bit  error  probability, 

Pi  =  p  ;  i  =  0,1... .,6-1 


(3.17) 

3  z“ 

Assume  that  for  each  word,  W,  a  fixed  amount  of  energy  ( Ew )  is  allocated. 
An  optimum  way  exists  of  asigning  the  energy  between  the  bits  which  minimizes 
the  m.s.e.  For  constant  amplitude  modulation  schemes,  the  way  of  assigning  dif¬ 
ferent  amounts  of  energy  between  the  bits  in  the  word  is  simply  by  having  different 
bit  duration  for  each  bit.  In  communication  systems  which  work  at  high  bit  rate,  this 
method  is  impractical  because  of  implementation  difficulties.  However,  for  systems 
with  low  bit  rate  (few  hundreds  of  bits/sec)  this  method  is  practical. 

In  Appendix  D,  the  optimal  bit  energy  asignment  is  calculated  for  a  DPSK 
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modulation  under  the  constraints  of  constant  word  energy.  Figure  3.11  shows  the 
optimum  normalized  bit  energy  Ebi  ;  i  =  which  gives  the  minimal  M.S.E 

as  a  function  of  the  signal-to-noise  ratio  for  a  12  bit  A/D  converter  and  a  binary 
DPSK  modulation  scheme.  The  normalized  bit  energy  is  defined  as  =  the  i-th  bit 
energy  normalized  by  E‘b  where  E‘b  is  the  total  word  energy  (Ew)  divided  by  number 
of  bits  (b). 

Normalized  Bit  Energy  For  12  Bits  A/D 


Figure  3.11  Bit  energy  alocation  that  gives  minimal  M.S.E  for  a  12  bits  A/D  converter 

and  BDPSK  modulation 


Figure  3.12  shows  the  ratio  between  the  M.S.E  when  all  bits  have  the  same  energy 
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and  the  case  when  energy  is  divided  optimally. 
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Figure  3.12  Comparison  of  the  M.S.E  when  equal  bit  energy  is  used  and  the 
optimal  M.S.E,  for  12  bit  A/D  converter. 


Another  intersting  way  to  achieve  unequal  bit  error  probability  is  by  using  a 
modified  error  correction  scheme.  This  class  of  error  corection  codes  (ECC)  family  is 
known  as  UEP  codes60.  They  have  the  property  that  a  code  be  designed  to  have  dif¬ 
ferent  error  correcting  capability  for  different  parts  of  the  word. 

In  their  pioneering  work,  Masnick  and  Wolf  60  have  introduced  the  concept 
of  unequal  error  protection  codes.  These  codes  have  the  property  that  some  of  the 
digits  in  a  codeword  will  be  decoded  correctly  if  /,  errors  occur,  while  others  will  be 
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decoded  correctly  only  if  f2  or  fewer  errors  ocurrs  where  /i>/2.  Boyarinov  and 
Katsman8  expanded  the  previous  work  and  described  a  broad  class  of  iterative  and 
concatenated  UEP  codes. 

IIL4  Do  we  need  error  correction? 

Inherently,  estimating  an  analog  waveform  from  its  digitized  samples  has  an 
estimation  error  caused  by  the  finite  step  size  of  the  A/D  converter  . 

This  error  usually  is  described  as  white  noise  and  is  known  as  quantization  noise67 
[Chapter  9]. 

When  digital  data  is  modulated  and  transmitted  through  a  non-ideal  channel, 
errors  occur  in  detecting  the  received  signals.  If  the  original  signal  is  an  analog  sig¬ 
nal  which  is  digitized  by  an  A/D  converter,  modulated  and  transmitted,  the  errors 
that  occur  during  the  detection  process  cause  errors  in  the  reconstructed  analog  sig¬ 
nal.  The  error  between  the  original  and  the  reconstructed  signals  is  a  noise-like  pro¬ 
cess  and  is  defined  as  the  detection  error  or  the  detection  noise. 

The  variance  of  the  A/D  quantization  error  is67: 


where  b  is  the  number  of  bits  of  the  A/D  converter.  For  the  detection  noise,  the 
mean  square  error,  e2,  (m.s.e)  is  given  by  (3.16)  (or  by  (3.17)  for  equal  bit  error  pro¬ 
bability).  If  we  keep  the  detection  mean  square  error  less  than  the  quantization 
error,  the  contribution  of  the  error  caused  by  wrong  decisions  of  the  receiver  will  be 
negligible  with  respect  to  the  quantization  noise  and  no  degradation  of  the  system 
performance  will  occur.  Let  us  require  that 
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2  -  2  • 
e2 


Substituting  (3.17)  and  (3.18)  into  (3.19)  leads  to  the  requirement  that 


12  2“(2“  -  1  )p. 


(3.19) 


(3.20 


P.< - =J - =  2~{2h+2)  . 

8  (2“  -  1) 

For  b  =  8,  has  to  be  less  than  1.9x10"*.  For  b  =  12,  pe  has  to  be  less  than  2.9xi0"8. 
For  optimal  bit  power  management  (III.3.3  and  appendix  C),  p,  must  be  less  than 
4xi0"5  for  b  =  8,  and  less  than  1.7x10"*  for  b  =  12 

Let  us  now  calculate  the  achieved  bit  error  probability  of  a  system  with 
source  level  of  170  dBi\\xPa,  where  the  distance  between  the  receiver  and  the 
transmitter  is  5  km 

Using  the  sonar  equation85,  the  received  signal  level  is 


EL  =SL-TL  +DI 


(3.21) 


where  the  transmission  loss,  TL,  is  given  by 


TL  =Ca  \og(R)  +  aR  . 

Ca  =  20  for  volume  spreading  and  Ca  - 10  for  cylindrical  spreading.  For  transmitter 
and  receiver  located  at  a  depth  of  few  hundred  meter,  ocean  depth  of  l  km,  and 
R  =5  km,  Ca  can  be  assumed  to  be  15.  For  a  =  1  dB/km ,  DI  =  1,  R  =  5  km  the  received 
signal  level  is  109  dB i\\xPa.  From  Urick85  [Figure  7.7,  curve  C],  the  sea  noise  level  at 
10  kHz  is  50  dB  l\\xPahT7z  . 

The  received  bit  energy  (for  bit  rate  of  1  kbit! sec)  is 

Eb  =  EL  T  =  109  -  30  =  79  [dBl\\iPa] 


%  =  =  29  dB 

No 


and 
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which  for  binary  DPSK  gives  a  bit  error  rate  much  lower  than  2.9xi0-8.  In  the  worst 
case,  when  Ca  =20,  the  transmission  loss  is  91  dB  which  gives  bit  error  probability  of 

l. 7x10^.  For  bit  a  rate  of  1  kbit / sec  and  an  8  bit  A/D,  the  achieved  bit  error  rate  (even 
for  the  worst  case  of  Ca  =  20  and  no  optimal  power  assignment)  is  on  the  order  of  the 
quantization  noise.  Therefore,  no  error  correction  (which  reduces  the  effective  bit 
rate)  is  required. 

m. 5  Combined  chirp  DPSK  and  error  correction 

As  we  saw  in  I3I.4,  no  error  correction  is  needed  for  the  given  scenario  in 
order  to  achieve  the  desired  system  performance.  However,  if  for  any  reason  the  sig¬ 
nal  to  noise  ratio  decreases  (reducing  the  transmitting  power,  higher  ambient  noise 
than  is  expected,  making  the  distance  between  receiver  and  transmitter  larger),  error 
correction  must  be  added  to  the  system. 

A  necessary  condition  for  chirp  DPSK  modulation  is  to  keep  the  time- 
bandwidth  product  of  the  signal  on  the  order  of  10.  Hence,  in  order  to  add  an  error 
correctic  capability  in  one  hand,  and  not  to  reduce  the  information  rate  on  the  other 
hand,  a  higher  order  DPSK  modulation  must  be  used.  In  Appendix  B  we  enclose 
curves  of  several  common  error  correction  techniques  [Odenwalder66  ].  From  these 
curves  and  from  Figure  3.9,  we  can  conclude  that  a  significant  coding  gain  can  be 

achieved.  As  an  example,  using  chirp  QDPSK  with  code  rate  of  ~  keeps  the  infor¬ 
mation  rate  and  the  time-bandwidth  product  the  same  as  chirp  BDPSK  with  the  same 
symbol  rate.  However,  if  a  bit  error  rate  of  10-6  is  desired,  a  signal  to  noise  ratio  of 
11.18  dB  is  required  when  BDPSK  modulation  is  use,  and  approximately  8  dB  is 
required  when  a  QDPSK  with  extended  (24,12)  Golay  code  is  used  (from  Figures 
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B.l  and  3.9).  When  a  convolutional  code  with  rate  of  and  constraint  length  of  7 

is  used  even  better  performance  can  be  achieved  (Figure  B.4  and  3.9)  and  a  signal  to 
noise  ratio  of  approximately  6.5  dB  is  required  to  get  the  same  performance  as 
before. 

If  a  higher  information  rate  is  required  (for  example,  2  kbit/sec  instead  of  1 
kbit/sec),  a  16- ary  chirp  QPSK  with  a  rate  y  error  correcting  code  can  be  used. 

When  using  the  previous  convolutional  code,  a  signal  to  noise  ratio  of  around  13  dB 
is  required.  Another  way  to  achieve  the  same  bit  rate  is  just  to  chose  a  chirp  QDPSK 
with  the  previous  convolutional  error  correcting  code  but  doubling  the  symbol  rate. 
In  this  case,  around  6.5  dB  signal  to  noise  ratio  is  required.  The  penalty  that  we  pay 
is  decreasing  of  the  intersymbol  interference  and  reducing  the  immunity  to  fading. 
This  is  due  to  the  fact  that  the  spectrum  of  the  chirp  signal  is  not  flat  any  more  and 
has  higher  side  lobes  on  the  other  hand.  Figure  3.13  shows  the  spectrum  of  a  chirp 
signal  with  time-bandwidth  product  of  5. 
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Figure  3.13:  Specra  of  chirp  signals  with  time-bandwidth  product  of  5  and  10. 


IV  Synchronization 


IV.  I  Introduction 


Power-efficient  digital  receivers  generally  require  a  digital  clock  synchron¬ 
ized  to  the  received  bit  stream  to  control  the  integrate-and-dump  detection  filter  or  to 
control  the  timing  of  the  output  bit  stream.  The  utility  of  any  particular  signaling 
waveform  depends  rather  critically  on  the  design  of  the  synchronization  system  used 
to  generate  the  receiver  carrier  phase,  the  fundamental  clock  and  bit  reference  from 
the  received  signal.  When  coherent  detection  is  used,  perfect  knowledge  of  the  car¬ 
rier  frequency  and  phase  is  vital  to  the  receiver  in  order  to  be  able  to  detect  the 
transmitted  signal.  For  incoherent  modulation  schemes,  only  the  knowledge  of  the 
carrier  frequency  is  necessary  to  detect  the  transmitted  signal. 

DPSK  modulation  is  semi  coherent  scheme.  Perfect  knowledge  of  the  bit 
duration  and  precise  bit  synchronization  are  necessary  at  the  receiver.  Synchroniza¬ 
tion  can  be  divided  to  three  levels  of  hierarchy: 

(1)  carrier  recovery  and  phase  synchronization. 

(2)  timing  recovery  or  symbol  (bit)  synchronization. 

(3)  word  or  frame  synchronization. 

First,  assuming  that  a  carrier-type  system  is  involved,  there  is  a  problem  of 
carrier  and  phase  synchronization  concerning  the  generation  of  a  reference  carrier 
with  a  phase  closely  matching  that  of  the  data  signal.  This  reference  carrier  is  used  at 
the  receiver  to  perform  a  coherent  demodulation  operation.  Next  comes  the  problem 
of  synchronizing  the  receiver  clock  with  the  baseband  data-symbol  sequence.  This 
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commonly  is  called  symbol  (or  bit)  synchronization.  Depending  on  the  type  of  sys¬ 
tem  under  consideration,  the  problem  of  word,  or  frame  synchronization  will  be 
encountered  further  down  in  the  hierarchy. 

Basically,  there  are  two  approaches  to  providing  synchronization.  Either  a 
separate  channel  is  allocated  for  sending  the  synchronization  signal  or  the  synchroni¬ 
zation  information  appears  on  the  same  channel  as  the  data  and  is  derived  from  the 
data  bearing  signal  itself.  As  discussed  here,  bit  synchronisation  is  restricted  to  self¬ 
synchronization  techniques  that  extract  the  clock  time  directly  from  the  noisy 
received  signal. 

When  a  coherent  modulation  approach  is  used,  precise  knowledge  of  the 
phase  and  the  carrier  waveform  is  necessary.  For  an  incoherent  approach,  scheme  the 
phase  of  the  signal  is  not  information  bearing  and  only  a  rough  knowledge  of  the  car¬ 
rier  waveform  must  be  known.  Both  approaches  need  good  bit  synchronization  in 
order  to  achieve  the  maximum  performance  of  the  system. 

For  DPSK  modulation  scheme,  the  information  bearing  signal  is  the  phase 
difference  between  two  adjacent  symbols.  As  the  signal  reference  is  always  the  pre¬ 
vious  pulse,  no  need  of  carrier  or  phase  recovery  exists.  But,  as  the  information  bear¬ 
ing  signal  is  the  phase  difference  between  two  adjacent  signal,  very  accurate  bit  syn¬ 
chronization  is  required.  By  bit  synchronization,  we  mean  time  of  arrival  and  bit 
duration  estimation. 

The  purpose  of  timing  recovery  in  DPSK  systems  is  to  recover  a  clock  at  the 
symbol  rate  or  a  multiple  of  the  symbol  rate  from  the  modulated  waveform.  This 
clock  is  required  to  adjust  the  delay  in  the  receiver  (Figures  3.7  and  3.8)  and  to  sam¬ 
ple  the  received  signal  at  the  correct  time  instants. 

Practical  timing  recovery  circuits  cannot  perfectly  duplicate  the  clock  used 


60 


at  the  remote  transmitter.  The  most  basic  requirement  is  that  the  average  frequency 
of  the  derived  timing  waveform  must  be  exactly  equal  to  the  average  frequency  of 
the  transmitted  signal.  Although  the  average  frequency  of  the  derived  timing  must  be 
exact,  the  timing  signal  usually  has  phase  or  timing  jitter.  In  DPSK  modulation  this 
timing  jitter  is  a  fundamental  impairment  and  dramatically  reduces  the  receiver  per¬ 
formance.  It  is  necessary  to  know  not  only  how  often  to  sample  the  data  bearing  sig¬ 
nal  and  the  exact  bit  duration  but  also  where  to  sample  it.  The  choice  of  sampling 
instant  is  called  the  timing  phase. 

For  conventional  modulation  approaches  (the  carrier  waveform  is  a  fixed  fre¬ 
quency  signal),  the  frequency  recovery  (which  enables  the  estimation  of  bit  duration) 
and  the  time  of  arrival  techniques  are  throughly  analyzed  and  well  covered  and  docu¬ 
mented  in  the  literature.  Simon  &  Lindsey54  (Chapter.  9),  Holmes45,  Ziemer  & 
Peterson97  (Chapter.  6),  Lee  &  Messerschmitt52  (Part  IV)  and  Franks  34  described 
and  analyze  the  most  common  methods  of  timing,  carrier,  and  bit  synchronization. 

Although  a  lot  of  effort  is  directed  toward  producing  stable  synthesizers  and 
clocks.  In  most  communication  systems,  a  long  term  instability  always  exists  and  a 
frequency  and  clock  shift  between  the  transmitter  and  the  receiver  occurs.  Carrier 
and  clock  recovery  techniques  which  synchronize  the  carrier  frequency  and  the 
clocks  ot  the  transmitter  and  the  receiver,  are  developed  and  analyzed  in  several 
references7’34,52,82.  In  most  communication  systems,  the  bit  duration  is  directly 
derived  from  the  carrier  frequency  and  synchronization  to  the  carrier  frequency 
enables  the  synchronization  of  the  bit  duration. 

Since  the  carrier  recovery  block  is  not  necesserlv  needed  for  DPSK  demodu¬ 
lation  and  all  the  common  methods  of  carrier  recovery  are  not  applicable  when  the 
carrier  is  a  chirp  signal,  a  different  scheme  of  bit  synchronization  is  proposed  here. 
This  method  is  based  on  the  fact  that  the  envelope  of  the  matched  filter  output  is  a 
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periodic  function  with  a  period  equal  to  the  bit  duration 

In  this  chapter,  we  will  present  methods  for  bit  synchronization  and  time  of 
arrival  estimation  for  a  chirp  DPSK  modulation  scheme.  Moreover,  we  will  show 
that  the  bit  synchronizer  is  not  sensitive  to  the  fading  and  the  multipath  phenomena 
which  are  dominant  in  underwater  acoustic  communication. 

IV.2  The  optimum  bit  synchronizer 


In  the  following  we  present  a  method  to  find  the  bit  duration  and  the  timing 
phase  for  chirp  DPSK  modulation.  The  received  signal  z(t)  [Figure  4.1]  is  character¬ 
ized  by: 

z(0  =  y(f,e)  +  n(r)=  £  y(t,e)  +  n(t)  (4.1) 

;  =o 

where  a,  takes  on  the  values  +A  or  -A  with  probability  y,  Raaj  =  A2  8(i  -  j),  and 

T  T 

y(t)  is  a  time  limited  signal  between  [-—  <  t  <  —],  which  depends  on  the  modulation 
scheme  used. 

Wintz  &  Luecke94  and  Lindsey  &  Simon54  showed  that  for  known  bit  dura¬ 
tion  T  and  for  a  time  interval  of  K  T  seconds  long,  the  maximum  a  posteriori  (MAP) 
estimator  of  e  is  the  e  which  maximizes 


K 

A (z  ,e)  =  £  In  cosh 


i=0 


~  j  z(t)y(t-iT-z)dt 
N0  r,(e) 


(4.2) 


The  ln(cosh(  ))  nonlinearity  is  approximately  a  square  law  device  for  small  argu¬ 
ments  and  has  approximately  a  magnitude  function  for  large  arguments,  that  is. 
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In  cosh(x )  = 


Ui 

2 

„2 


UI  »  1 
UI  «  1 


(4.3) 


Using  the  approximation  in  (4.3),  a  suboptimum  bit  synchronizer  configuration, 
which  can  easily  be  implemented  is  given  in  Figure  4.1. 


Figure  4.1  Bit  synchronization  -  block  diagram  of  a  suboptimum  scheme. 

The  noisy  input  bit  stream,  z(t),  is  first  fed  to  a  filter  matched  to  the  input 
wave  shape  (y(r)X  then  squared  and  fed  to  the  bandpass  filter.  The  bandpass  filter 
output  is  approximately  a  sine  wave  whose  positive-going  zero  crossings  estimate 
the  data  transition  times.  The  bandpass  filter  bandwidth  determines  the  syncronizer 
memory  and  the  variance  of  the  timing  jitter.  A  limiter  and  an  antibouncing  circuit  is 
used  to  provide  a  rectangular  shaped  periodic  signal  which  drives  the  timing  circuits 
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(Figure  4.1). 

The  output  of  the  matched  filter,  z,(r),  is 

zi(0  =  yi(^)  +  "i(0  (4-4) 

where 

yi('  e)  =  iaiy1(t-iT-E)  (4.5) 

i=0 

and 

>i(0  =  y(0  *  MO 
/* i (r )  =  MO  *  h(t)  . 

h(t)  is  the  impulse  response  of  the  matched  filter.  After  squaring  and  low  pass  filter¬ 
ing,  (the  lowpass  cutoff  frequency  coq  <  co^p  <  2  co0 ),  we  get 

Z3(t,e)  =  y30  £)  +  n-i(‘)  ,  (4-6) 

where 

y3  (r  ,e)  =  A2  £  y3(f-tT  -  e)  +  DC  term 

i  =o 

and  n3(t)  is  a  zero  mean  random  process.  y3(r  ,e)  also  can  be  written  as 

y3  (t  ,e)  =  A2 y3(t)  *  £  8(t-iT -e)  +  DC  term  . 

;=o 

The  power  spectrum  of  y3  (0  is 

y3(co)  =  A  2exp(-;  C0£)  f  3  (CO)  X  0>r  2  5(a>  -  i  a>r)  +  DC  term  (4.7) 

i=0 

where  cor  =  -y- 

Since  the  spectrum  of  y3(co)  consists  only  of  discrete  frequencies  which  are  multiples 
of  the  bit  rate,  bandpass  filtering  y3(co)  with  a  narrow  filter  centered  at  one  of  the 
discrete  frequencies  gives  the  timing  phase  and  the  bit  rate  or  any  multiple  of  the  bit 
frequency.  If  the  bandpass  filter  is  centered  around  cor,  then 


y*(t)-AT  cos(coT  +<?>) 


(4.8) 
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where 


Aj  —  (Or  A2  Y 3(0)7 ) 

and  y3(o)  is  the  Fourier  transform  of  y3(/). 


IV.2.1  System  performance 

Two  metrices  generally  are  used  to  measure  the  performance  of  a  bit  synchronizer, 
'he  rms  jittei  or  timing  misalignment  and  the  mean  of  the  absolute  value  of  the  jitter. 
In  this  chapter,  we  will  calculate  the  rms  jitter  of  the  proposed  scheme.  No  general 
expression  of  timing  misalignment  exists  and  diffemt  results  exist  for  any  bit  syn¬ 
chronization  implementation. 

Wintz  &  Luecke  94  developed  results  for  a  RC  lowpass  filter  with  3-dB 
banwidth  B  =  y  followed  by  a  square-law  detector  used  to  syncronize  a  NRZ  PAM 

modulation  scheme.  The  clock  component  is  recovered  with  a  phase  locked  loop  or  a 
bandpass  filter  with  an  equivalent  bandpass  memory  of  KT  (KT  is  the  observation 
time).  The  expected  magnitude  of  error  for  a  raised-cosine  input  wave  form  is 
approximately  equal  to 

l  0  33 

^£[|e|]  =  -^r  ;  Yj,  >  5,  K  >  18 

T  SK  yb 

L.  E.  Franks34  ,  Meyers  &  Franks61  ,  Holmes45  and  Chiu  &  Lee20  also  have 
analyzed  different  types  of  bit  synchronizers  and  presented  their  performance. 

The  signal  component  at  the  matched  filter  output  is  y\(t)  =  y(t)*h{t)  and  the 
noise  componnent  n{(t)  is:  «,(/)  =  where  h(t)  is  the  matched  filter  impulse 

response  and  is  equal  to:  h(t)  =  y'(-t).  The  (•)*  denotes  the  complex  conjugate  of  (■). 
Without  losing  generality,  e  can  be  assumed  to  be  zero.  Then  for  y(t)  a  narrow  band 
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process,  the  match  filter  output,  (r  t(r ))  can  be  written  in  quadrature  form  as 


zi(f)=|  ^aigit -iT)  +  nc(t)  |  cos(cooO  -  (0  sin(o)of )  ,  (4.9) 

i=° 

where  g(t)  is  the  lowpass  equivalent  match  filter  output  to  the  input  y(t). 

The  lowpass  filter  output ,  z3(0>  is  then 

1  4  r  1  'l 

(4.10) 

zi(t)  is  a  nonstarionary  noncentral  chi-squared  random  process  with  two  degrees  of 
freedom  and  probability  density  function  ,  p2,(z3,D,  equal  to 


1 

2 

z3  (t)='A 

I,  digit  ~iT)  +  nc(t) 

+  !/2  ns(t) 

L 1=0 

L 

Pz,(z  3.0  = 


2a2. 


exp 


[(J3(f))2  +  Z3]2 


2a. 


Vzi (0 


Chd)))2 


a2. 


U(zz) 


and  (  •)  =  (  )  modulo  T.  The  autocorrelation  R.y(t,t-x)  of  z3(r)  is: 

R„(U-x)  =  E  [z3(0z3(i-t)1 


(4.11) 


(4.12) 


A  ^  — 


g2(Og2(r-x)  +— j-'  J  °'-2, 


2g 2(t  —  t)  +  2g2(7)  f  +  A  2/?„,( t)  g  (f)g  (T-t)  +  a„,  -  2/?  3  (x) 


Define 


where 


«3(0  =  *3(0- >3(0 


y3(f)  =  £  [z3(r)]  =  y/t2g2(f)  +  CT2  . 


(4.13) 


For  n(0,  a  AWGN  process  with  two  sided  spectral  density  —  and  .4  as  defined  in 
(4.1),  nz(i )  is  a  zero  mean  random  process  with  autocorrelation  Rn,(t,x)  given  by 

r  t2 


-N0 

Rn,(t,x)  =  E  [  ni(t)nj(t  -  x)  }=  A2—g(x)  g(T}g (T-x)  + 


Nr 


-g(  t) 


(4.14) 


For  a  large  signal  to  noise  ratio,  n3(t)  is  well  approximated  by  a  non  station¬ 
ary  zero  mean  Gaussian  random  process  with  autocorrelation  function  as  in  (4.14), 
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and 


cr2 n,(t)  =  A2  g2(T)<j1n  +  a4„,  .  (4.15) 

n4(0  is  also  a  zero  mean  Gaussian  random  process  with  autocorrelation  function  is 

equal  to 


R*lt,D=Rnlt,x)*  hB At)*  h * bj> (~t ) 
where  hBP(t)  is  the  impulse  response  of  the  bandpass  filter. 

The  estimated  timing  phase,  e,  is  the  instant  when  zA{t)  crosses  the  zero  axis 
(Figure  4.2). 


* 


Figure  4.2  Bit  synchronizer  -  estimation  error. 

Define  u  as  the  error  between  the  true  e  and  the  estimator  (e).  Without  loss  of 
generality,  assume  that  e  =  0.  The  variance  of  the  bit  synchronizer  error  is  then: 

T 

2 

o2  =  J  u2  P  (o)  d u  . 

JL 

2 
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Assuming  that  the  autocorrelation  function  of  n4(t)  is  differentiable  and  using  the 
analysis  of  the  well-known  level  crossing  problem  (Papoulis68  Chapter  11-4)  we  get 
(assuming  an  antibouncing  circuit  exists  in  the  limiter  block) 


T 


2  1 

1  c  2 

'  -«"».(0,0)  1 

1 

2 

[Ar  sin  (  &>/•  -u)  ]2 

—  J  v2 

71  Jr 

*-.(0,0) 

exp 

2  7?n4(0,0) 

2 


(4.16) 


='  1  ”2 

TT  J 


j  m25„4(0,m)  4  co 

J  SJO, co)4co 


exp 


[At  sin  (  o>r  u)  f 

2  j  S..< 0.0.)  |f  J 


4v> 


where 


*"J0.0)=-J^  KJf.i) 


r=fl,-c=0 


For  an  ideal  bandpass  filter  with  Bbj,  bandwidth  and  a  channel  free  from  multipath, 
substituting  (4.8)  and  (4.15)  in  (4.16)  we  get 


2  _±_ 


J  o^S^O.to)  4(0 
J  5J0, (0)4(0 


(4.17) 


x  J  u2exp 


o>r  Yfc2  5^  (ay)  sin2(o>r  u) 


2  Y* 


2 


8k2[-£-*(0)  +  1]  |  5f‘(co)  4 


(0 


2 


4u 


4  2  T 

where  5 .•(©)  is  the  Fourier  transform  of  g2(0,  and  y*  is  defined  as  — —.  The  value 

'  2Nn 


of 


I  (O^.to.o))  4(0 
J  S„4(0,(o)  4(0 


is  known  as  the  r.m.s.  bandwidth  of  the  process  n4(0  and 
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denoted  by  co2l  2  (Helstrom41  pp  474). 

Since  the  bandpass  filter  is  relatively  narrowband,  we  can  assumme  that  in  the 
relevant  bandwidth  Sx>(co)  is  constant  and  is  equal  to  o>r) .  In  this  case, 

B,, 

2 

J  Sgiu)d(ti  =  BB_p  SgK(Or)  (4.18) 

Bar 

2 

and  for  ideal  banpass  filter 


J  co2S„4(0,co)  d  to 


J  S„4(0,co)  d  co 


For  a  chirp  input  signal  defined  as: 


y  ( t  )  =Re  [  exp  j  ( co0 1  +  lA|u2 )  ]  y<i  <  ^ 

The  matched  filter  impulse  response  h  (t)  is 


h  (r )  =  Re  [  exp  j  ( co<,  t  -  'A\U2 )  ]  -y  <  t  <  y 

The  signal  componnent,  y?l  (r),  at  the  matched  filter  output  is 


y\ (0  =  COS(COoC  )  g  (t) 

where 


2(0=  yy  sin(y  [fj-T  |  r  |  -  r2])  ;  -T  <t  <T 


(4.19) 

(4.20) 

(4.21) 
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Time  Normalized  to  Bit  Duration 


Figure  4  J  The  equivalent  lowpass  matched  filter  output  -  g(t) 


Using  "  K  "  as  defined  by  Wintz  &  Luecke  94,  Bbj>  can  be  aproximated  as 


K 


After  some  simple  substitutions  we  get 


2_JL 


co 


j  u2  exp 


K  coj  7*  S^^cor)  sin2(cor  u) 

8  7T2  [  +  1  ] 


dx>  . 


(4.22) 


Figure.  4.4  gives  the  normalized  error  standard  deviation  -y  as  a  function  of  the  sig¬ 
nal  to  noise  ratio  for  different  values  of  K  when  the  chirp  signal  has  a  time- 
bandwidth  product  of  20n  and  bit  rate  (y)  of  1  kbits/sec 
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S l gpo  1  fo  Noise  Ratio  (  DB  ) 

Figure  4.4  The  normalized  error  standard  deviation  as  a  function  of  the  signal  to  noise  ratio 

for  different  values  of  K. 

IV.III  Bit  synchronization  using  time  of  arrival  estimation  method 

In  the  ocean  acoustic  channel,  there  are  two  major  obstatcles  towards  the 
successful  use  of  the  bit  synchronizer  scheme  presented  above.  The  frequent  non- 
specular  multipath  implies  noise  in  addition  to  the  ambient  noise,  and  the  presence  of 
specular  multipath  does  not  guarantee  tracking  of  the  principal  timing  phase.  In  the 
presense  of  multipath,  the  signal  component  at  the  matched  filter  output  is 

yi('.e)  =  ya(^)  +  yw(»,e)  +  y„(rx)  ,  (4.23) 

where  yd(tx)  is  the  direct  path  component  and  is  equal  to 

yd(t  .£)  =  Z  a,  yi(f-iT-e) 

i=0 

y ™(r  ,e)  is  the  specular  multipath  contribution  and  is  given  by 
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y«(<  .£)=  L  ^  yi(t-iT-E-  0 

i=0 

where  C  is  the  time  delay  between  the  direct  and  the  indirect  paths.  For  a  stationary 
or  a  very  slowly  moving  receiver  and/or  transmitter,  the  time  delay  £  is  a  slowly 
varying  function  with  a  time  varying  mean  and  a  small  variance. 

The  nonspecular  multipath  contribution  y„(/,£)  is 
y«  (f  .£)  =  2  c-  >  i _jT_e  +  9« ) 

i= 0 

where  0,  is  uniformly  distributed  between  0  and  2n.  For  this  case,  the  output  of  the 
square  law  device  (Figure  4.1)  is 

y3  (0  =  IA2U‘  -e)  +  k2U‘  ~  01  *  £  8  (t-iT)  +  C2y3(0  *  £  8 (r-iT  +  0,)  (4.24) 

i=0  i=0 


+  21 
i= 0 


a,y3(t  -  iT  -£)&,> 3(r-(T-Qcos(a-£)  +  a,y3(t-iT-£)c,y3(r-jT-0,)cos(£-0,) 


+  bty3(t-iT-Qc,y3  (r-jT-^cos^-©,)  j 

where  a  is  the  phase  difference  between  the  direct  and  the  indirect  paths,  A  =  |  a,  | 
and  B  =  \bi\  =  kA  for  all  /.  The  first  term  in  (4.24)  is  a  periodic  function  with  period 
T.  The  second  and  the  third  terms  are  noise-like  signals.  In  the  frequency  domain 


y3(co)=A2 


exp  (-j  uz)  +  k2  exp  (-j  o)  a) 


Y3  (u)  x  ajj-  £  8(co  -  i  a>r)  +  NOISE 

i=0 


Figure  4.5  Bit  synchronization  -  spectrum  of  the  low  pass  filter  output 
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It  can  be  seen  that  multipath  is  a  destructive  phenomenon  for  the  bit  syn¬ 
chronization  concept  described  above.  It  increases  the  overall  noise  level  and 
prevents  good  bit  synchronization.  The  estimated  timing  phase  e  falls  somewhere 
between  the  real  e  and  a  which  is  far  away  from  e. 

Modification  of  single  pulse  time  of  arrival  estimation  is  an  alternative  bit 
synchronization  concept  which  enables  us  to  synchronize  to  the  timing  phase  even  in 
the  presence  of  specular  multipath.  As  before,  assume  that  the  bit  duration  T  is  per¬ 
fectly  known.  In  such  a  case  it  is  easy  to  prove  [  Helstrom41  page  376  ]  that  the  max¬ 
imum  likelihood  estimate  e  of  e  is  the  time  at  which  the  rectified  output  of  the 
matched  filter  is  maximum.  The  output  will  have  many  peaks  and  the  highest  of  them 
identifies  the  time  e.  Figure  4.6  shows  a  block  diagram  of  the  system  for  the  multiple 
case. 


Z4O)  =  y  4(0  +  «4(0  2  3(0  =  y  3(0  +  « 3(0 


Figure  4.6  Bit  synchronization  -  autocorrelation  implementation 


The  synchronizer  output,  z4(t,e),  is 


24«,e)  =  -^£z3(r-zT) 

"  0 

In  the  presence  of  multipath,  the  signal  component  at  the  syncronizer  output  is 


(4.25) 


The  third  term  in  (4.26)  is  due  to  the  nonspecular  multipath  phenomenon  and  is  com¬ 
pletely  random  in  phase  and  amplitude.  The  forth  term  is  also  a  random  process.  For 
0,  uniformly  distributed  between  -n  and  n,  these  two  terms  have  a  zero  mean  and  for 
K  large  enough,  they  vanish. 

For  a  single  pulse,  the  variance  of  the  estimator  e  is41 

Var  i=  [  2  y^P2]"*  (4.27) 

and  for  multiple  ulses, 

Vare=[  2Ky4p2 

where 


j  [j;y(t)]2]dt 
B2  =  — - 

r  -H- 

|  [y(t))2)dt 

or,  according  to  Parseval’s  theorem 


J  [o>2  |  Y( co)  \2]d(o 

p2=- - 

r  -H* 

J  [  |  Y( co)  |  ]2]4co 

For  incoherent  detection  (which  is  our  case)  y(t),  is  substituted  by  g{t). 


1  K~l 

The  signal  z4(0  =  —  ^  z3(r  -pT)  and  explicitly, 

K  p=o 


Ut)=±A282(T)+±Vc  A)2(t-Q  +  ~ 


Z  + 


Z  n2(t-pT) 

P- 0 


(4.28) 


1  AT-l  - 

+  7  I  I  g[t-(i+p)r]  g[f-(i+p)7'-Qcosa,_p 

A  p=0  .=o 
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+  g[*~(i+P)T] 


+  TtS  ncV-pT)  £  *,-pg  [/ -</  +p  )T -Qcosa,_p 

^  n 


+  TrZ  n>(t-pT)  X *.-p5[Mj+p)r-^]sina1_p 

K  P  i 


Fory4(/)  defined  as;y4(D  =  £[24(D]  and  n4(t)  defined  as  n4(/)=24(/)-y4(/)  we  get 


y4(t )  =  —  A2g2(D  +  —  k2A2g\t-Q  +  +  kA2g (Dg (f-Q£ [cosa]+a2  .  (4.29) 


and 


Rn.(t  J-t)  =  E  [« tU  )n  M  -r)]  = 


(4.30) 


77  XL  R»:  -p  >r +*> +  4r  -p  >r +T> 

p  ?  pi 


g  (P)g  (i -~(r +Q)E  [cosap  ] 


+  k  g  (t  -t)  g  (t -£+(£-t ))£  [cosa,,  J  +fc2  g  (/  -Q  $  (r  -<-t)£  [cosap  cosa,  ] 


+  k  g  (t ~Q  g  </-£-r)£ [sinap sina?  ]  +*2g  (Dg (T-t)  +  o‘*  A2 g(t)g  (T-Q  £ fcosap  ] 


+  FA4g(7)g(7-i;)g(r  -^-t)  g  U -C-r)[  £  [cosa] 
where  r( ,  =/>  -  #  of  bits  in  ( • )  . 


l] 2  -  < 


For  /?„,  =  -y-  g(7),  (4.30)  becomes: 


tF= 


(4.31) 


1 

No 

2  ,  AT0  -  ' 

~K 

ygh) 

+  42  —  g(T) 

r.^.  — 

g  (D  g  (T-t)  +  ■  -  k  g(T)g  (i  ~(j+C))E  [cosa] 


r—  _  _  _  _ 

+  —jjr-  k  g(»-0g(^+(^-7))£lcosa]  +k2  g(t-Q  g(t-^-r) 


75 


C  ' 

+  k2A* g(D  g(7-C)  gd-^+t^-T)  g(t-C,-T)  cosa  2 

*  « 

In  the  second  case,  the  phase  a  varies  very  rapidly  and  is  asummed  to  be  uniformly 
distributed  between  Ht-Jt].  Then: 

>4(')=yAV(7)  +  |  k2A2g2(t- 0  +~g(0)  (4.34) 

and 

=  -~giT)  +A2-y-g(T)[  g(Dg(T-T)  +  k2gU-Qg(r^-7)  (4.35) 

In  the  presense  of  specular  multipath,  the  performance  of  the  estimator  is  not 
only  governed  by  the  additve  noise,  but  also  by  the  probability  of  locking  onto  the 
indirect  path  instead  of  the  direct  one.  Let  us  define  the  case  when  the  syncronizer  is 
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locked  to  the  peak  at  e  as  State  1,  and  when  it  is  locked  on  the  peak  at  t,  as  State  2. 
The  locking  probabilities  state  diagram  is  depicted  in  Figure  4.7 

P  i/i  P  m 


Figure  4.7  Error  state  diagram. 


If  the  bit  synchronizer  is  in  State  1 ,  it  will  not  change  its  state  in  the  next  observation 
if  the  peak  at  e  is  higher  than  at  C,.  The  probability  for  this  event,  P  m  ,  is 

+oo  x 

P i/i/oc  =  J  J  PUQUJ®UM=xj^)=q]dqdx  (4.36) 

Xm—ao  jo 

For  large  signal  to  noise  ratio  and  large  K,  Pu 0U4O  can  be  assumed  as  normal  and 
then 


P[n/a  = 


°Mo>  onJQ  Vl  -p2 


(4.37) 


+«>  x+Ay 

x  |  |  exp 


<?  +  Qhort 1 

2  [  Om0)O^q  l2  (  1  -  p2) 


dq  dx 


where  Ay  =  y4(0)  -y4(Q  and 


P2  = 


«,.( 0.0 


°-140)°n40 

After  some  simple  algebraic  manipulation,  we  get 


Pm  = 


o^o) 

_ -fao 

2  J  exp 

\  ^ 

4> 

.  2  o240) 

L 

Ay[onJ0)-po^0] 


dr  . 


(4.38) 


Assume  the  synchronizer  was  initially  locked  onto  the  direct  path  (State  1). 
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If  a  decision  is  made  every  T,  seconds,  (TS=K  T  where  K  is  the  number  of  bit  aver¬ 
aged),  the  mean  time  of  being  locked  on  State  1  is 

Tlxkl  =  TsPm  P'n  .  (4.39) 

vv 

The  mean  time  of  being  locked  onto  state  2,  if  initially  locked  on  State  1  is 


1  lock  2  ' 


T'PmP 


2/1  1,2  (^7 


Since  Pv2  =  1  -  Pm,  and  P 2/1  =  1  -Pn  we  get  that 


(4.40) 


and 


Ttock  1  ~ 


PmTs 

'~Pm 


P lock!  ~ 


TsPvi 
'-P  2/2 


(4.41) 


A  = 


Tiocki  __  P  \n('-P 2/2) 
Tlockl  =  (l-/>,/,)2 


(4.42) 


Figures  4.8-4. 10  shows  the  dependence  of  log  A  on  the  signal  to  noise  ratio  (yb)  for  K 


=  100  and  for  several  values  of  k  and  a,  for  the  case  when  a  varies  very  fast. 


Figure  4.8  A  as  a  function  of  signal-to-noise  ratio  for  various  values  of  k,  and  for  a=C. 


LOG  LRMBDfl 
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Stgnol  io  Noise  Rcrtio  C  DB  ]  Stgnol  io  Noise  Rortlo  [  OB  ] 


0  i  8  12  16  20 


Signal  to  Noise  Ratio  [  DB  ] 


Figure  4.9  A  as  a  function  of  signal-to-noise  ratio  for  various  values  of  k,  and  for 
«--0.25rc. 
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Figure  4.10  A  as  a  function  of  signal-to-noise  ratio  for  various  values  of  k,  and  for  a=n. 


Figure  4.11  shows  the  dependence  of  log  A  on  the  signal-to-noise  ratio  (yb),  for  K  = 
100  and  for  several  values  of  k  for  the  case  where  a  changes  very  rapidly. 


LOG  LAMBDA 


V  Equalization 


V.l  Introduction 

The  most  severe  problem  in  underwater  acoustic  communication  is  inter¬ 
symbol  interference  (ISI).  The  major  causes  of  the  intersymbol  interference  are:  a) 
channel  multipath  and  b)  shape  of  the  modulated  data.  A  careful  inspection  of  the 
character  of  the  underwater  communication  channel  shows  that  two  kinds  of  ISI  exist 
-  that  due  to  "macro"  multipaths  and  that  due  to  the  "micro"  multipaths.  Delayed 
replicas  of  the  original  signal  (caused  by  reflection  from  the  ocean  surface  and  bot¬ 
tom,  which  can  be  on  the  order  of  several  thousand  of  bits)  are  scrambled  at  the 
receiver  with  the  signal  received  from  the  direct  path  and  are  disastrous  to  the 
receiver.  The  "micro"  multipaths  result  in  spreading  the  signal  in  addition  to  produc¬ 
ing  ISI  from  some  of  the  adjacent  symbols. 

Even  though  multipath  seriously  degrades  the  receiver  performance,  it  is 
known  that,  the  theoretical  channel  capacity  of  a  wide  band  link  is  not  seriously 
reduced46.  This  theoretical  result  explains  why  so  much  effort  is  directed  toward 
solving  the  channel  equalization  problem.  However,  this  result  applies  only  when 
exact  knowledge  of  the  channel  parameters  is  exit  which  is  a  condition  that  hardly 
ever  is  met  in  practice.  Therefore,  adaptive  receiver  structures  are  required  in  order 
to  suppress  interference  arising  from  unknown  and/or  changing  multipath  channels. 
Considerable  research  has  been  performed  investigating  performance  criteria  for 
optimizing  the  equalizer  parameters.  Since  the  most  meaningful  performance  meas¬ 
ure  for  a  digital  communication  system  is  the  average  probability  of  error,  it  is  desir¬ 
able  to  choose  the  filter  parameters  to  minimize  this  performance  index.  However, 
the  probability  of  error  is  a  highly  nonlinear  function  of  the  equalizer  parameters. 
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Consequently,  the  probability  of  error  as  a  performance  index  for  optimizing  the  the 
equalizer  parameters  is  impractical. 

Two  criteria  have  found  widespread  use  in  optimizing  equalizer  parameters. 
One  is  the  peak  distortion  criterion  and  the  other  is  the  least-mean-square  error  cri¬ 
terion.  The  peak  distortion  criterion  is  defined  as  the  worst-case  intersymbol 
interference  at  the  output  of  the  equalizer.  The  minimization  of  this  performance 
index  is  called  the  peak  distortion  criterion.  The  least-mean-square  error  (LMS)  cri¬ 
terion  is  defined  as  minimizing  the  mean- square  value  of  error  between  the  equalizer 
output  and  the  desired  signal.  The  LMS  equalizer  minimizes  the  signal  to  distortion 
ratio  at  its  output  within  the  constraints  of  the  equalizer  time  span  and  the  delay 
through  the  equalizer. 

In  general,  the  equalizer  is  a  recursive  filter  with  poles  as  well  as  zeros.  A 
filter  with  poles  has  two  serious  disadvantages  compared  with  an  all-zero  filter:  (1)  it 
becomes  unstable  if  the  poles  move  outside  the  unit  circle  (during  the  adaptive  pro¬ 
cess,  or  because  of  round  off  errors)  and  (2)  its  performance  surface  generaly  is  non- 
quadratic  and  may  have  local  minima92.  For  these  reasons,  recursive  adaptive  filters 
are  found  only  in  a  limited  number  of  applications.  These  disadvantages  and  the  fact 
that  any  transfer  function  can  be  expressed  as  an  all-zero  transfer  function  with  an 
infinite  number  of  zeros  are  the  reasons  why  most  equalizers  are  implemented  as 
transversal  filters  (where  filters  with  an  infinite  number  of  zeros  are  approximated  by 
filters  with  a  finite  number  of  zeros).  Both  RF  as  well  as  undersea  channels  are  time 
varying.  Therefore,  adaptive  equalizer  structures  are  required  in  order  to  suppress 
interference  arising  from  unknown  and/or  changing  multipaths. 

Much  work  has  been  done  in  the  field  of  equalizer  performance  and  imple¬ 
mentation  both  as  fixed  as  well  as  adaptive  equalizers.  The  paper  by  Qureshi76  and 
books  by  Proakis75  (Chap.  6)  and  Lee  and  Messerschmitt52  (Chap. 8)  provide  a  good 
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review  of  the  equalization  literature.  The  books  by  Haykin39  and  Widrow92  provide 
good  review  of  adaptive  filtering  concepts. 

However,  none  of  the  existing  equalizer  configurations  can  be  implemented 
directly  for  an  underwater  acoustic  channel.  Two  characteristics  of  the  underwater 
acoustic  channel  makes  the  existing  equalizer  structures  impractical:  (1)  the  very 
long  multipath  on  the  order  of  some  thousands  of  bits  and  (2)  the  non-minimal  phase 
character  of  the  channel.  These  two  basic  characteristics  of  the  underwater  acoustic 
channel  makes  it  completely  different  from  a  conventional  radio  link.  The  reason 
why  the  delay  of  the  multipath  is  so  large  is  due  to  the  speed  of  sound  in  the  ocean 
being  relatively  small.  The  nonminimal  phase  characteristics  of  the  channel  is  due  to 
the  multipath  from  the  sea  surface  is  often  being  much  stronger  than  the  direct  path. 

When  a  zero  forcing  criterion  is  used,  the  transfer  function  of  the  equalizer  is 
given  by75 

(51) 

where  //c(z,r)  is  the  z  transform  of  the  channel  impulse  response.  When  the  LMS  cri¬ 
terion  is  used,  the  transfer  function  of  the  equalizer  satisfies75 


(5.2) 


Hc(z,t)+N0 

In  both  cases,  the  transfer  function  of  the  equalizer  is  proportional  to  the  inverse  of 
the  transfer  function  of  the  channel.  Due  to  the  large  delay  of  the  multipath  arrivals, 
the  transfer  function  of  the  channel  has  many  zeros  and,  as  a  result,  the  transfer  func¬ 
tion  of  the  equalizer  has  many  poles.  As  mentioned  previously,  the  equalizer  must  be 
implemented  as  an  all  zero  transfer  function  due  to  stability  problems.  Implementing 
the  equalizer  as  an  all  zero  transfer  function  requires  a  number  of  taps  which  is  much 
larger  than  the  maximal  delay  of  the  channel.  When  a  short  multipath  exits,  this  con¬ 
dition  is  reasonable,  but  when  the  delay  is  very  long  as  in  the  case  of  ihe  underwater 
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channel,  such  an  implementation  is  not  acceptable.  Furthermore,  a  more  serious 
problem  is  the  non-minimal  phase  property  of  the  channel  since  a  stable  inverse  of  a 
non-minimal  phase  transfer  function  does  not  exist. 

V.2  An  LMS  equalizer  for  underwater  acoustic  channel 


Following  (2.9)  the  z  transform  of  the  channel  impulse  response  is: 

He  (Z,t)  =  Oo  +  £  Ct;  Z~D‘  (5.3) 

1=1 

where  Dt  is  the  number  of  samples  at  7D;.  Without  loss  of  generality,  it  can  be 
assumed  that  the  largest  a  is  ak  k  *  0.  In  this  case,  (5.3)  can  be  written  as  : 

Hc(z,t)  =  oikrD‘HXz,t)  (5.4) 

where 


m-k 


Hc(z  ,/)=  £  °i  z 


'  -D, 


t=>-k 


,  a, 

a,  =  —  and  D;  =D,-Dk.  Furthermore  we  can  assume  that  ak  =  1 .  In  this  case,  the 
a* 

strongest  path  is  considered  to  be  the  main  path  and  obviously,  H'(z,t)  is  a  noncausal 
filter.  Instead  of  writing  H'(zj)  as  a  summation,  it  can  be  written  as  a  product: 

N  L  N-L 

#c(z)  =  Ho  n  (Zi  -z~1)  =  H0n  ( z!-z~l)  n  (5.5) 

i=0  1=0  1=0 

z\  are  the  zeros  of  Hjz )  which  are  inside  the  unit  circle  and  z°  are  the  zeros  which  are 
outside  the  unit  circle.  The  index  t  has  been  omitted  for  convenience  but  it  is  clear 
that  the  transfer  function  is  time  dependent.  For  high  signal  to  noise  ratio,  the 
transfer  function  of  the  equalizer  reduces  to  be  the  inverse  of  the  channel  transfer 
function  no  matter  what  optimization  criterion  is  used: 
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H,  (z )  =  -7—  =  H‘{z )  H?(z ) .  (5.6) 

Hc(2) 

Hg(z )  and  H?(z)  are  the  causal  and  the  noncausal  parts  of  Ht(z),  respectively, 
where 


and 


Hct{2)  = 


1 

Hfr) 


(5.7) 


l 

H?(z) 


(5.8) 


H?(z)  is  a  noncausal  filter  and  cannot  be  implemented  practically.  Rewriting  (5.8)  as 


//~(z)  =  S  3"  z’"  =  zA  (TA  X  p"  z-  (5.9) 

and  approximating  the  infinite  sum  by  a  finite  sum  of  A  terms,  delaying  by  z~A  and 
substituting  into  (5.6)  leads  to 


Ht(z)=— - l- - 2  3"  z-"  (5.10) 

n  (z/-z_1)  ns:0 

i=0 

He{z),  as  appears  in  (5.10)  is  a  stable  causal  filter  and  can  be  implemented. 

Due  to  the  problems  of  instability  and  convergence,  it  it  preferable  to  imple¬ 
ment  a  filter  (especially  if  it  is  adaptive)  by  a  FIR  filter.  Any  HR  filter  can  be  approx¬ 
imated  by  a  transversal  filter  with  a  finite  number  of  taps.  A  rule  of  thumb  says  that 
the  required  number  of  taps  is  5  to  10  times  the  number  of  poles  in  the  IIR  filter. 
When  the  number  of  poles  is  relatively  small,  such  implementation  is  acceptable.  On 
the  other  hand,  when  the  number  of  poles  is  very  large  (typically  a  few  thousand  in 
the  underwater  communication  channel),  implementing  a  transversal  filter  with  the 
number  of  taps  5  to  10  times  larger  is  not  practical. 

An  alternative  method  to  directly  inverting  the  channel  transfer  function 
must  be  found.  Since  almost  perfect  knowledge  of  the  original  data  exits  at  the  out¬ 
put  of  the  decision  device,  we  are  enable  to  build  the  equalizer  as  a  modified 
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combination  of  a  feed  forward  equalizer  and  a  version  of  a  decision  feedback  equal¬ 
izer  52,75  (DFE).  The  idea  behind  this  configuration  is  to  use  the  strongest  path  as  the 
reference  rather  than  just  the  first  arrival.  The  received  signal  y(n)  is  sent  to  the 
feed-forward  filter  h?{n)  (Figure  5.1)  which  filters  out  the  noncausal  part  of  the  input. 
The  output  x(n)  is  the  causal  part  of  y(n).  Now,  the  data  stream  taken  from  the  out¬ 
put  of  the  decision  device  d(n)  is  sent  to  the  feedback  filter  hct{n)  which  approximates 
the  causal  part  of  the  channel.  Subtracting  this  output  u(n)  from  x(n)  will  leave  us 
with  the  strongest  path  signal.  The  architecture  in  Figure  5.1  looks  like  more  like  an 
echo  canceller  than  a  classic  equalizer.  The  strongest  signal  is  almost  always  the 
main  path  or  the  first  multipath.  Thus,  the  noncausal  filter  in  (5.8)  has  relatively  few 
poles  and  can  be  realized  as  a  transversal  filter  (the  summation  pan  in  (5.10)). 


Figure  5.1  A  block  diagram  of  an  underwater  channel  equalizer. 

If  d(n)  is  the  signal  transmitted  by  the  transmitter,  the  received  signal  y\n)  is 


y\n)  =  d(n)  *  hc(n) 


(5.11) 


88 


where  hc(n)  is  the  impulse  response  of  the  channel.  Defining  y(n)  =  y\n  +Dk),  taking 
the  z  transform,  and  substituting  (5.5) 


Y (z )  =  D  (2 )  He (2 )  =  D  (z )  Hcc(z )  H ~(z ) 


(5.12) 


From  Figure  5.1: 


X(z)=Y(z)H?(z) 

where  (•)  denotes  the  estimator  of  (•).  If  H?(z)=  — - —  ,  we  get  that: 

//*(z) 


X(z)=D(z)HI(z). 

On  the  other  hand,  the  output  of  the  decision  device  is  fin)  where: 

/(n)  =  £4  5(n)  (5.13) 

and  d,  is  the  detected  signal  and  takes  on  the  values  +1  or  -1.  For  all  practical  pur¬ 
poses,  it  can  be  assumed  that  d,  =dt  for  all  i.  I(n)  is  sent  to  the  filter  hT(n)  which 
generates  the  signal  d{n)  .  d(n)  is  the  detected  version  of  the  original  transmitted  sig¬ 
nal,  and  in  high  signal  to  noise  ratio,  d{n)~d(n).  The  sequence  d(n)  is  fed  into  h^(n), 
which  is  the  feedback  part  of  the  equalizer  where: 

f  Kin )  n  >  i 

*«">-{  0  else  *514> 

The  z  transform  of  the  output  of  Kin )  C/(z)  is: 

Uiz)  =  Diz)  Hc,iz)  =  Diz)[Hcciz)  -  h‘(0) )  (5.15)) 

=  Xiz)-hW)Diz) 

where  Diz)  is  the  z  transform  of  din).  Tne  z  transform  of  the  output  of  the  equalizer 
D(z)  is: 

Diz)  =  Xiz)-Uiz)  (5.16) 

=  Diz)h‘iO) 


Note  that  if  ^(0)  =  1,  then  din)  =  din) 
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V.2.1  r  uiculation  of  filter  weights 

In  developing  the  filter  weight  calculation  method,  it  is  convenient  to  use 
matrix  notation.  The  composition  of  the  Mcxl  tap- weight  vector  of  hce(n)  is  defined  as 

iv- : 

wj  =  [wf  ,  W%, .  VVw.  ] 

The  tap-weight  vector  of  hct(n)  is  defined  as  : 

w^  =  [wf  . Wffm] 

The  superscript  T  signifies  transposition.  The  A^xl  output  vector  at  time  n  is  defined 

as  d(n): 

dT(n)=[d(n),d(n- 1) . d  (n  -Mc  + 1 )] 

and  the  M^xi  input  vector  at  time  n  is  defined  as  y(n): 

yT(n)=[y(n),y(n-l) . yi.n-M^.+l)] 

Define  the  estimation  error  e(n)  as: 

e(n)  sd(n)  -  d(n)  (5.17) 

where  d(n )  is  the  detected  value  and  d(n)  is  the  estimated  value  of  d(n)  at  the  input 
to  the  decision  device.  The  mean  squared  error  is: 

J=E(\e(n)\2)  (5.18) 

where  £(•)  is  the  expected  value  of  (■).  Writing  (5.18)  explicitly  one  gets: 

J  =e\  [  vv"  y(n)-w”  d(n)-d(n)  ]  [yH  (n)  wv  -  dH{n)wc  -  d'  (n)  ]  (5.19) 

After  some  simple  algebraic  manipulations  we  get: 

7  =  vy^  ^  vv^:  -  w’J  Rdy  wK  -  w. 


(5.20) 


90 


(5.21) 


(5.22) 


PT  ^\pjy,-pj] 

Rj 

R  = 

~{L*y  Rd 

Minimizing  the  mean  square  error  leads  to  the  well  known  Wicner-Hopf  equation: 

W„=R-'P  (5.23) 

When  the  tap-weight  vector  equals  its  optimum  value  W0,  the  mean  square  error 
attains  its  minimun  value  Jmin  which  is  equal  to 

Jnun^^-P"  R']  P  (5.24) 

One  approach  to  solving  (5.23)  is  by  using  direct  matrix  inversion  tech¬ 
niques.  Although,  this  procedure  is  quite  straightforward,  it  can  present  serious  com¬ 
putational  difficulties  (especially  when  the  filter  contains  a  large  number  of  tap 
weights  and  the  input  data  rate  is  high  as  is  in  our  situation).  An  alternative  pro¬ 
cedure  is  to  use  the  method  of  the  steepest  descent63  which  leads  to  the  following 


-PdyWnc-  Pjy  +  'ALc  Hd  +  Pd  +  w"  p"  +  CT 

where  the  superscript  H  signifies  Hermitian  transposition  and 

Ry  =  E\y_(n)  y" (i>)], 

Rdy  =  E[d(n)yH(n+A)], 

R_d  =  E[d(n)  dH (n)], 

Pd=E  [d{n)  dH (n)], 
pJy=E{d(n)f(n+A)] 

Rearranging  (5.20)  leads  to  the  familiar  form: 

J=WhRW-PW-PhWh+G2 

where 

WT  *[w£,weT] 
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recursive  relation: 


V^(n  +  l)=iy(/i)  +  4[F -R  W(»)],  n  =0,1,2,...,  (5.25) 

The  parameter  n  controls  the  size  of  the  incremental  correction  applied  to  the  tap- 

weight  vector  as  processed  from  one  iteration  cycle  to  the  next. 

Usually,  the  values  of  R  and  P  are  unknown  and  have  to  be  estimated.  The 
most  common  way  of  estimating  R  and  P  is  to  use  instantaneous  estimates  that  are 
based  on  sample  values  of  the  tap-input  vector  and  desired  response  as  defined  by: 


and 


*00  ■ 


y(n)f(n) 
-d(n)yH  (n  +A) 


-y  ( n +A) dH ( n ) 
d(n)dH(n) 


(5.26) 


PT(n)  =  [d(n)yH(n),  d(n)dH(n)]  (5.27) 

Substituting  (5.26)  and  (5.27)  in  (5.25)  will  lead  to  the  well  known  LMS  adaptive 
algorithm  proposed  by  Widrow  and  Hoff91. 


V.3  Error  propagation 


In  the  decision-feedback  equalizer  (DFE),  the  equalizer  depends  on  results 
produced  by  the  decision  device  in  the  receiver.  When  no  error  exists  in  detecting  the 
transmitted  data,  the  DFE  always  has  an  advantage  (and  often  a  substantial  advan¬ 
tage)  over  purely  linear  equalization31.  However,  when  perfect  detection  of  the 
transmitted  signal  does  not  occur,  the  presence  of  error  propagation  due  to  the  effect 
of  past  decision  errors  raises  the  question  of  whether  error  propagation  is  catas¬ 
trophic  and  if  not,  over  what  interval  does  the  effect  last  on  the  average. 

Following  the  definitions  of  Figure  5.1, 
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d(n)  =  x(n)  -  u(n)  (5.28) 

where 

x{n)  =  d{n)*  h£(n)  +  nc(n)=I(n)  *  hT(n)*  h‘{n)  +  nc{n)  (5.29) 

nc(n)  is  the  ambient  noise  passed  through  the  filter  h™(n),  and: 

u(n)  =  d(n-Ns )  *  fan)  =  l\n-Ns)  *  hT{n)  *  fan)  (5.30) 

Using  the  channel  model  defined  by  (2.9),  substituting  in  (5.13)  and  (5.14),  and 

TDi 

using  the  assumption  that  — —  is  an  integer  number  of  bits  for  all  i ,  one  obtains: 

*  S 

m— £ 

d(n)  =  dH  hT(n)+  £  a,  idn-TD.  ~  ^h-td,)  hr(n)  (5.31) 

■=o 

TD 

The  assumption  that  equals  an  integer  number  of  bits  for  all  i  corresponds  to  the 

*  S 

case  that  maximal  damage  occurs  from  a  wrong  decision.  The  DFE  equalizer  is 
applicable  only  in  systems  which  have  relatively  small  errors  at  the  decision  device 
output.  Therefore,  one  can  assume  that  only  one  error  in  the  m-k  relevant  bits 
occurs.  The  worst  case  is  when  an  error  exists  in  the  dt  which  corresponds  to  the  a  of 
the  largest  absolute  value.  Without  loss  of  generality  the  error  can  be  assumed  the 

n  -  1'*  bit. 

Eight  different  cases  can  occur  as  indicated  in  Table  5.1. 
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Case  No. 

dn 

dn- 1 

4„-i 

1 

1 

1 

1 

2 

1 

-1 

1 

3 

1 

-1 

-1 

. 4~ 

1 

1 

-1 

5 

-1 

1 

1 

6 

-1 

-1 

1 

7 

-1 

-1 

-1 

8 

-1 

1 

-1 

Table  5.1:  Possible  combinations  of  present,  previous,  and  detected  bits. 


For  symmetry  reasons,  Cases  1-4  are  the  same  as  Cases  5-8. 

For  Cases  1  and  3  : 

d  (n)  =  d  (n)  +  nc 

and  the  corresponding  bit  error  probability  is 

Eb 

P.l=P.3  =  f^) 

No 

where  the  function  /(•)  depends  on  the  modulation/demodulation  scheme  is  used. 
For  Case  2: 


and 


d(n)  =  d{n)  +  2  cti  d(n)  +  nc 


Pt2=f[(l  +  2a,)2 -f] 
Nn 


and  for  Case  4: 
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d  (n )  =  d  (ft )  -  2  ot]  d  (n )  +  nc 

with 

P<4=/[(l-2ai)2^] 

yv0 

Again,  due  to  symmetry,  the  same  results  are  accepted  for  a  positive  or  a  negative. 
Without  loss  of  generality,  let  us  assume  that  a  >  0. 

The  relation  between  making  an  error  or  correct  decision,  depending  on  the 
previous  decision  is  a  two  state  Markoff  process  as  depicted  in  Figure  5.2. 


Figure  5.2  The  relationship  between  the  decision  of  the  present  and  the  previous  bits 

Given  that  the  previous  bit  was  detected  incorrectly,  then  the  present  bit  will  be 
detected  incorrectly  with  probability  P,2  for  the  events  described  as  Cases  2  and  8  in 
Table  5.1.  The  present  bit  will  be  detected  incorrectly  with  probability  PtA  for  the 
events  described  as  Cases  4  and  6  in  Table  5.1. 

The  average  time  between  an  error  decision  and  the  next  correct  decision  Ttp 
is: 
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Ttp  =  |  [2  -  ( P,2  +  Pt4 )]  £  n  [|  ( P<2  +  P,4  )]»  - 1  (5.32) 

_ _ 1 _ 

1- |  (^2  +  ^4) 

In  the  worst  case,  Pt 4  =  1  and: 

J -  (5.33) 

Since  fep  is  a  finite  number,  the  error  propagation  is  of  finite  duration.  Moreover, 
since  Pel  is  a  very  small  number,  the  average  error  propagation  is  no  more  than  two 
bits.  Hence,  under  condition  of  low  bit  error  rate,  the  degradation  of  system  perfor¬ 
mance  due  to  the  error  propagation  is  negligible. 

V.4  Numerical  example 


Assume  a  channel  with  the  following  impulse  response: 


hc  {n )  =  d  (n )  -  2.6d  (n  - 1)  +1 ,25d  ( n  -2)  -  0.  Id  (n  -3)  (5.34) 

and  transfer  function  Hc(z) 

Hc(z)  =  1  -  2.6z-1  +  1.25z-2  +  O.lz-3  =  ( 1  -  2z"!  )(1  -  0.5Z'1 )(  1  -  (Hz'1 )  (5.35) 

Hc{z)  has  three  zeros,  two  of  them  are  inside  the  unit  circle  and  one  (z  =  2)  is  outside 
the  unit  circle,  which  causes  Hc(z)  to  be  an  non-minimum  phase  transfer  function. 
Following  the  definitions  in  (5.7)  and  (5.8),  we  obtain: 


///(z)  = 


1 


1 


"?(:)■ 


Hcc(z) 

1 


(1  -  0.5z~‘  Xl-O.lz"1) 


(5.36) 


1 


-1 


-1-m 


//^(z)  (1  -  2z-1 ) 


3TT  =  I  [2z-‘] 


Choosing  A  =  5  and  following  (5.9)  and  (5.10): 
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«.(*) 


_ 1 _ 

(1  -  0.5z-1  Xl-O.lz-1) 


2-5  z5  £  [2z-1]" 

n=0 


From  (5.8)  and  (5.14)  we  obtain: 


and 


2~5  £  [2z-'r 


n=0 


H‘(z)  =  -0.6+0.05z-1 


(5.37) 


(5.38) 


(5.39) 


A  computer  simulation  was  run  to  illustrate  the  performance  of  the  method. 
A  long  sequence  of  +1  or  -1  randomly  selected  was  sent  to  a  block  with  transfer 
function  as  in  (5.34)  which  simulated  the  under  water  acoustic  channel.  The  output 
y  (n )  of  this  block  was  sent  to  a  block  which  simulated  the  receiver  as  in  Figure  5. 1 
and  the  weights  of  the  two  filters  were  calculated.  In  order  to  illustrate  the  perfor¬ 
mance  of  the  method,  a  single  pulse  was  transmitted  and  analyzed  in  the  receiver 
with  the  filter  weights  that  had  been  calculated  before.  The  example  of  the  single 
pulse  shows  how  the  equalizer  cancels  all  the  multipaths  except  the  strongest  one. 
Figures  5.3  -  5.6  provide  simulation  results  where  a  single  impulse  was  transmitted 
by  the  transmitter,  passed  through  a  channel  with  the  impulse  response  in  (5.34)  and 
processed  by  the  receiver  in  Figure  5.1.  In  each  figure,  the  weights  of  the  filters  were 
calculated  in  a  different  way.  Figure  5.3  corresponds  to  filter  weights  calculated 
directly  from  (5.38)  and  (5.39). 
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y(n) 


x(n) 


u(n) 


din) 


5AMPLE5  SAMPLES 


Figure  5  J  Equalizer  weights  calculated  directly  from  (5.38)  and  (5.39) 

yin)  is  the  distored  signal  received  by  the  receiver,  x(n)  is  the  output  of  the  forward 
filter  hj(n),  u(n )  is  the  output  of  the  feedback  filter  h‘(n),  and  d(n)  is  the  equalized 
signal.  Figure  5.4  shows  the  results  for  the  case  where  the  filter  weights  were  calcu¬ 
lated  using  the  Wiener-Hopf  equation.  In  both  cases,  the  channel  was  considered 
noiseless.  Figure  5.5  and  Figure  5.6  show  the  results  when  the  LMS  algorithm  was 
used.  In  Figure  5.5,  the  channel  was  assumed  to  be  noiseless  and  in  Figure  5.6,  a 
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noisy  channel  (AWGN)  was  considered  where  —  =  10  dB. 

W0 


y(n) 


u(n) 


x(n) 


d(") 


Figure  5.4  Equalizer  weights  calculated  using  the  Wiener-Hopf 
LMS  algorithm  in  noisless  channel. 


RKfLJTUDE  AMPLITUDE 
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y(n) 


SABLES 


x(n) 


u(n) 


5fl^LE5 


d(n) 


5AVPLE5 


Figure  5 .5  Equalizer  weights  calculated  using  the  LMS  algorithm 
(500  iterations)  (no  noise). 


AMPLITUDE  AMPLITUDE 


VI  Examples  and  Discussion  of  Results 


A  new  technique  to  combat  destructive  channel  fading  for  underwater  acous¬ 
tic  communication  has  been  proposed.  This  technique  combines  both  chirp  and  dif¬ 
ferential  phase  shift  keying  (DPSK)  modulation.  Chirp  DPSK  modulation  is  an 
efficient  modulation  scheme  in  the  presence  of  a  slowly  time-varying,  frequency 
selective  fading  channel. 

Ignoring  the  fading  problem  and  neglecting  its  destructive  effects  causes  a 
communication  link  which  behaves  nicely  in  a  non-fading  environment  to  be  useless 
under  fading  conditions.  The  bit  error  probability  which  is  proportional  to  the  inverse 
exponent  of  the  signal-to-noise  ratio  in  a  non-fading  environment,  becomes  inversely 
proportional  to  the  signal-to-noise  ratio  itself  under  fading  conditions.  Thus, 
dramatic  degradation  in  system  performance  occures.  As  an  example,  achieving  a 
bit  error  probability  of  KT6  in  binary  coherent  PSK  with  no  fading  requires  a  signal- 
to-noise  ratio  of  approximately  10.5  dB.  In  the  presence  of  a  Rayleigh  fading  chan¬ 
nel,  achieving  the  same  bit  error  probability  a  requires  54  dB  of  signal-to-noise  ratio. 
Using  frequency  diversity  methods  improves  the  performance  of  the  communication 
link.  The  more  parallel  branches  that  are  used,  the  better  the  system  performance. 
The  random  nature  of  the  fading  phenomena  makes  a  coherent  modulation  scheme 
useless  since  the  channel  itself  destroys  the  coherence  of  the  transmitted  signal. 

In  the  presence  of  slowly  time-varying  (relative  to  the  bit  rate),  frequency 
selective  fading,  the  chirp  DPSK  modulation  is  a  good  anti-fading  modulation 
scheme.  Chirp  modulation  is  very  attractive  as  an  anti-fading  modulation  scheme  due 
to  its  spectral  properties.  For  large  time-bandwidth  products  ( BTe ),  the  power  spec¬ 
trum  has  almost  a  rectangular  shape,  (i.e.  a  flat  spectrum  in  the  relevant  bandwidth 
and  almost  zero  outside  of  it).  The  flatness  of  the  spectrum  is  similar  to  using  a  large 
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number  of  parallel  diversity  branches.  The  equivalent  number  of  branches  is  on  the 
order  of  the  number  of  degrees  of  freedom  (K)  of  the  channel.  The  maximum 
number  of  equivalent  diversity  branches  is  achieved  without  increasing  the  problem 
of  cross-talk.  Another  advantage  of  chirp  DPSK  modulation  is  the  simplicity  of  the 
receiver  realization  which  is  equivalent  to  a  single  channel  DPSK  receiver. 

When  dealing  with  frequency  diversity  techniques,  the  available  number  of 
diversity  branches  depends  on  the  ratio  between  the  available  bandwidth  B  and  the 
bit  rate  Rc  (the  bit  rate  is  the  inverse  of  the  bit  length  Tc).  Assuming  the  information 
bandwidth  is  approximately  Rc,  then  in  order  to  prevent  strong  cross-talk  between 

adjacent  channels,  the  number  of  diversity  channels  is  bounded  above  by  L  <  —  . 

On  the  other  hand,  the  number  of  degrees  of  freedom  (K)  in  the  chirp  DPSK  modula¬ 
tion  technique  depends  only  on  the  nature  of  the  channel  (i.e.  on  the  ratio  between 

D 

the  available  bandwidth  B  and  the  frequency  selectivity  pH  where  K  =  — ).  The  chirp 

P  H 

modulation  technique  can  be  considered  plausible  only  for  cases  when  the  time- 
bandwidth  product  (BTC)  is  greater  than  or  equal  to  10.  Under  this  condition,  the 
assumption  of  a  flat  spectrum  and  constant  phase  of  the  chirp  signal  across  the 
relevant  bandwidth  holds.  The  spectrum  and  the  phase  of  a  chirp  signal  with  BT  =  10 
is  shown  in  Figure  6.1. 
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Figure  6.1  The  spectrum  and  the  phase  of  a  chirp  signal  with  BT  =  10 


As  an  example  of  the  performance  of  the  chirp  DPSK  modulation  scheme 
relative  to  other  frequency  diversity  methods,  let  us  consider  a  shallow  water,  short 
range  underwater  acoustic  communication  system.  The  available  bandwidth  is  10 
kHz  and  the  required  data  rate  is  1  Kbit/sec.  From  measured  data18,  the  frequency 
selectivity  p*  is  approximately  50  Hz.  For  this  example,  the  time-bandwidth  product 
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(BTC )  is  10000  Hz  x  1  msec  =  10  ,  L  <  10,  and  K  =  =  200. 

In  Figure  6.2  we  present  a  comparison  of  various  modulation  schemes  with 
10  diversity  branches  and  a  chirp  DPSK  with  20  and  200  degrees  of  freedom.  For  all 
practical  purposes,  200  degrees  of  freedom  can  be  considered  as  infinite. 
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Figure  6.2  Comparison  between  various  diversity  methods  and  chirp  DPSK 

modulation 

Observing  the  piots,  we  see  that  in  the  presence  of  slowly  time-varying,  fre¬ 
quency  selective  fading,  using  chirp  DPSK  is  superior  to  all  other  modulation 
schemes  using  frequency  diversity  and  a  gain  of  almost  10  dB  in  signal  to  noise  ratio 
is  achieved  with  respect  to  NBFSK  with  diversity  of  10. 

Another  way  to  combat  the  fading  channel  is  to  use  combinations  of  coding 
and  modulation  schemes.  This  method  optimaly  combines  diversity  modulation  and 
coding.  A  modified  version  of  the  combined  coding  and  modulating  approach  that 
was  proposed  by  Pieper  at  el71  is  used  by  a  group  at  Woods  Hole  17  in  their 
telemetry  system. 
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Comparing  the  performance  of  the  chirp  QPSK  modulation  with  bandwidth  expan¬ 
sion  of  5  with  several  other  combined  coding  and  modulating  schemes  (Figure  3.10) 
having  similar  bandwidth  expansion  (4),  shows  that  the  QDPSK  is  superior  over  the 
other  methods.  A  necesery  condition  for  any  DPSK  scheme  to  be  used,  is  that  a  high 
coherency  between  two  adjacent  symboles  must  exist.  In  the  case  of  slowly  time- 
varying  fading  and  quasi-static  terminals,  the  coherency  between  two  adajecent  sym¬ 
boles  is  almost  one  and  the  degradation  in  the  system  performance  is  negligible  (Fig¬ 
ure  A.  8) 

An  underwater  channel  for  acoustic  communication  u  a  very  tough  channel 
in  the  sense  of  having  disastrous  interference  caused  by  strong  multipath  signals.  A 
necessery  condition  for  reliable  communication  is  to  find  methods  which  eliminate 
the  effect  of  the  multipaths.  Although  the  results  of  sea  experiments  show  that  only 
one  dominant  path  exists,  it  is  a  mistake  to  assume  that  this  is  the  case  in  general. 
Hence,  a  way  to  prevent  the  huge  ISI  caused  by  the  multipaths  has  to  be  found. 

Multipath  has  two  major  properties  that  makes  the  equalizer  implementation 
very  difficult: 

(1)  The  delay  between  the  direct  path  and  the  multipath  signal  can  be  on 
the  order  several  thousand  bits  and  hence  an  equalizer  with  many 
taps  is  required. 

(2)  The  intensity  of  one  of  the  multipath  signals  can  be  much  higher 
than  the  direct  path  signal,  i.e.  the  channel  is  a  nonminimum  phase 
channel. 

In  this  thesis,  we  propose  an  equalizer  implementation  which  requires  a  minimum 
number  of  taps  and  overcomes  the  problem  of  the  nonminimun  phase  of  the  channel 
impulse  response.  The  proposed  equalizer  is  a  combination  of  forward  and  decision 
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feedback  equalizers.  In  this  configuration,  the  strongest  signal  (which  is  not  neces¬ 
sarily  the  direct  path)  is  assumed  to  be  the  reference  signal  and  hence  the  channel 
changes  to  being  a  noncausal  instead  of  a  nonminimum  phase  channel.  The  forward 
filter  approximates  the  part  of  the  equalizer  that  equalizes  the  noncausal  part  of  the 
channel  and  the  feedback  filter  equalizes  the  causal  part  of  the  channel.  The  decision 
feedback  configuration  enables  implementation  the  equalizer  with  the  minimum 
number  of  taps.  Although  a  nonconventional  equalizer  configuration  is  chosen,  it  is 
shown  that  using  the  mean  square  error  optimization  criterion  leads  to  the  well 
known  Wiener-Hopf  equation  for  the  tap  weights.  An  error  propagation  phenomena 
exists  but  it  proves  to  be  non-catastrophic  with  an  average  error  propagation  time  of 
2  bits.  Some  numerical  simulation  results  are  included  which  demonstrate  the  perfor¬ 
mance  of  the  proposed  equalization  method. 

The  strong  multipath  phenomena  results  in  bit  synchronization  being  com¬ 
plicated.  In  order  to  achieve  the  best  performance,  the  synchronization  method  was 
designed  in  such  a  way  that  the  synchronizer  is  always  locked  to  the  strongest  mul¬ 
tipath.  The  intensity  of  the  multipath  signals  change  with  time  (especially  the  micro 
multipath)  so  in  order  to  prevent  the  synchronizer  from  losing  lock  and  locking  onto 
another  multipath,  an  average  multipath  intensity  is  calculated  and  the  synchronizer 
is  locked  onto  the  maximum  average  multipath  signal. 

The  telemetry  system  proposed  and  analyzed  in  this  work  is  not  a  universal 
system.  It  is  important  to  emphasis  that  the  system  performes  optimally  only  under 
specific  condition  (source/reciever  almost  at  rest  and  a  slowly  time  varying  channel). 
When  the  channel  is  not  a  slowly  time-varying  channel,  the  propopsed  methods  of 
equalization  and  bit  synchronization  do  not  work  any  longer  and  different  methods 
must  be  found.  Designing  a  telemety  systems  for  a  fast  time  varying  channel  is  a 
tough  problem  which  has  been  not  solved  yet.  Since  the  number  of  the  mutipath 
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signals,  their  intensity  and  location  varies  rapily,  no  reference  exists  for  synchroniza¬ 
tion.  Moreover,  the  rate  of  the  multipath  variation  is  so  fast  that  the  common  adap¬ 
tive  equalization  algorithm  do  not  succeed  in  following  the  variation  of  the  channel. 
New  adaptive  algorithms  with  a  fast  adaptation  properties  must  be  found.  If  the 
requirement  of  general  telemety  system  exists,  the  system  probably  will  be  much 
more  complicated  (and  expensive)  and  will  have  worse  performance  than  a  system 
that  has  been  designed  to  operate  in  a  specific  well-defined  scenario. 

VI.  1  Real  data  analysis 

A  stream  of  pseudo  randomly  selected  ones  and  zeros  was  differentially 
modulated  by  a  chirp  signal  and  transmitted  as  a  part  of  the  at-sea  experiments.  A 
receiver  implemented  in  software  processed  the  data.  Figure  6.3  shows  a  block 
diagram  of  the  receiver. 


V, 


Bandpass 

Filter 


Whitening 

Filter 


■e 


Y  =S  +JV, 


Match 

Filter 


Bit 

Synchronizer] 


1  bit 
Delay 


K5>- 


Integrate 
and  Sample 


Figure  6.3  Implemented  DPSK  receiver  -  block  diagram. 
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The  ocean  ambient  noise  is  not  white  and  has  a  power  spectrum  shown  in 
Figure  6.4. 


Figure  6.4  Power  spectral  density  of  the  ambient  noise  (after  bandpassing) 

It  can  be  approximated  roughly  as  a  single  pole  AR  process.  Therefore,  the  whiten¬ 
ing  filter  was  implemnted  as  a  simple  derivative.  The  received  signal-to-noise  ratio 
(yb)  was  quite  high  (20  dB).  In  order  to  analyze  the  performance  of  the  system,  at 
several  values  of  signal-to-noise  ratio  ambient  noise  (N2  measured  at  different  time 
interval  which  was  clean  from  any  signal)  was  added  artificially  to  the  received  sig¬ 
nal  (y). 

A  stream  of  6xi06  bits  was  analaysed.  The  achieved  bit  error  and  the  theoret¬ 
ical  bit  error  rate  for  DPSK  modulation  in  the  presence  of  AWGN  are  presented  in 
Figure  6.5. 
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Figure  6.5  Comparison  between  theoretical  and  achieved  bit  error  rate. 

As  shown  in  Figures  2.4,  2.5  and  2.6,  the  first  multipath  is  much  stronger  than  the 
direct  signal  and  the  other  multipath  signals.  Therefore,  the  bit  synchronizer  was 
locked  onto  the  first  multipath.  Since  the  first  multipath  was  so  large  relative  to  the 
other  signals,  in  the  range  of  bit  error  rate  that  was  analyzed,  the  absence  of  the 
equalizer  did  not  affect  the  achieved  results. 


Appendix  A  -  Chirp  as  a  Frequency  Diversity  Technique 


Introduction 

The  main  idea  behind  all  frequency  diversity  methods  is  to  spread  the  spec¬ 
trum  of  the  signal  as  much  as  possible  over  the  transmission  bandwidth.123  This  is 
precisely  the  reason  why  spread  spectrum  communication  is  a  useful  diversity  tech¬ 
nique. 

A  simple  and  very  useful  spreading  method  is  chirp  carrier  signal  modula¬ 
tion.  As  with  all  familiar  frequency  diversity  techniques,  it  is  strongly  desired  to 
make  use  of  as  many  branches  as  possible  for  a  given  frequency  bandwidth.4 
Although  the  channel  and  the  branches  are  designed  (at  least  theoretically)  to  be 
orthogonal,  choosing  the  branches  too  close  together  in  the  frequency  domain  results 
in  of  cross-talk  and  intersymbol  interference  between  adjacent  branches. 

A  chirp  signal  is  the  same  as  an  isolated  signal.  Thus,  no  cross-talk  between 
adjacent  branches  occurs.  Moreover,  the  rectangular  character  of  the  chirp’s  power 
spectrum56  (at  least  for  time-bandwidth  products  much  greater  than  one)  minimizes 
the  intersymbol  interference  problem.  In  this  appendix  we  will  analyze  the  perfor¬ 
mance  of  FSK  and  DPSK  modulation  schemes  using  a  chirp  as  the  carrier. 

A.l  Possible  signal  sets  for  binary  chirp  modulation 

In  binary  chirp  modulation,  when  both  mark  and  space  occupy  the  same  fre¬ 
quency  region,  three  signal  sets  can  be  chosen.  These  are  labeled  (a),(b),  and  (c)  and 
are  given  by 

(a) 
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Ill 


5,(0  = 


sin'1 


cos 


(OW  + 1/2& ) 


(A.l) 


5°(0  =  ((Hot  -  Vi\xx2)  . 


Here,  the  correlation  coefficient,  p,  between  the  two  signals  S,(0  and  50(O  is 
defined  as: 


vi  r 


J  5,(0  S0(t)dt 


P=T 


-ViT 


—  if 

TT  9  L 


j  Sxit)  dt  j  si(t)dt 


2  -'AT  L 


cos  ( 2oty  )  cos  (  p/2) 


dt  .  (A. 2) 


(b) 


Sl(f)=[co"](tDo,-'^2) 
5°(0=  (sfn)  _  ^MJ2)  • 


(A. 3) 


The  correlation  coefficient,  p ,  between  S,(0 ,  and  S0(O  is  : 


nr 


dt  . 


(A.  4) 


P  =  -r  J  I  cos  ( 2o)o/  )  sin  ( pi2 ) 

Z  0  L 

For  appropriate  values  of  p  and  coq  ,  p  vanishes  and  the  two  signals  are  orthogonal, 

(c) 


5,(0  =  [cos]  (“o'  +  ‘Ap/2) 
5°(0  =  [sf„]  (“o'  -  ‘/ip/2)  . 


(A. 3) 


In  this  case  : 


p  =  ^y  J  sin  (p/2 )  dt  . 


(A. 6) 
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A.2  Detection  of  incoherent  BFSK  chirp  modulation  in  the  presence  of 
AWGN  and  Rayleigh  fading 

Before  analyzing  the  performance  of  binary  chirp  modulation  in  the  pres¬ 
ence  of  AWGN  over  a  fading  channel,  let  us  start  with  the  simple  case  of  a  channel 

distorted  only  by  additive  white  Gaussian  noise  with  two  sided  spectral  density  . 

For  the  following  signal  set,  the  receiver  is  a  incoherent  BFSK  receiver  as 
shown  in  Figure  A.  1 . 


Figure  A.l  Incoherent  BFSK  Receiver 
In  this  figure, 

S\(t)  =  A  co^coo t  +  Vip/2),  ~y  <  t  <  —  (A.7) 

and 

S0(t)  =  A  cos(coo t  -  Vip/2),  -y  "  1  -  \  ■ 

Without  any  loss  of  generality,  it  can  be  assumed  that  St(t)  is  transmitted.  It  is  easy  to 
show  (neglecting  the  high  frequency  terms)  that 
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and 


Then 


*4>-f 


r  0.5  T 

U0(-z)  =  A  J  cos(p  t2)dt 
1  0 


2  NoT 
°N  =  —r- 


(A.8) 
(A. 9) 


and 


P(ji/5l)  = 


y\ 

exp  • 

y?  +  uhj) 

•  /0- 

yxUt  y) 

on 

2  o% 

On 

(A.  10) 


yo 

P(yo!si)  =  -r  exp 
On 


yo  +t/o(y) 


2a£ 


/n 


y0u0(j) 


On 


Define  the  probability  of  correct  decision  as  Pe  .  Then  [Helstrom41  p.  431] 


(A.ll) 


where: 


PC=  J  J  P(y\/Si)P(yo/Si)dyidy0 

y«=0  yt=yt 

=  G|Vyi,V%|  -  Vi  exp  [ -Vi  ( Yo  +  Yi)l  2o  I'ffOo] 


(A.  12) 


VH  f) 

7‘  "  2^ 


i  =0,1 


and  Q  (a,P)  is  the  Marcum  Q  function.  The  bit  error  probability  Pe  is  given  by 


Pt  =  \-Pc  ■ 


114 


A.3  The  Incoherent  case  in  the  presence  of  frequency  selective, 
slowly  time- varying  Rayleigh  fading 


Due  to  the  nature  of  the  fading  channel,  the  transmitted  signal  is  strongly 
distorted  in  amplitude  and  phase  when  it  arrives  at  the  receiver.  This  distortion 
renders  a  coherent  system  absolutely  useless  in  a  strongly  fading  environment.  For 
this  reason,  we  present  only  the  analysis  of  a  incoherent  system. 

For  the  signal  set  in  (A.7),  the  received  signal  is 


y  i(t )  =  Re  y  i(t )  exp  0  GW )  r  . 


(A.  13) 


where  (see  (1.3)) 


y1(r)  =  aAS1(0+  J  ASj(/)  h(t-x:t)  dz+n(t) 


(A.  14) 


Sx(t)  =  exp  (Vi/ \u2)  ;  — <  t  < 


The  signal  component,  U\(~) ,  at  the  output  of  the  receivei  is: 


Ul{2)=2  J  \  )$!(')+  J  ASMhlt-rjiidxSUt)  \  dt  .  (A.15) 


where  ( •  )*  denotes  the  complex  cojugate.  For  convenience,  define  r(t)  as: 


r(t)=  A  j  Sx(x)  h(t-v,tx)  dz  . 


r(t)  also  can  be  expressed  as: 


(A.  16) 


r(t)  =  A  j  Sj(w) //(©;*!)  exp  O' ©0 


(A.  17) 


=  A  J  1 S  i (to)  1  exp(j  8(w))  H  (co;/  0  exp  (/  co/ ) 
For  a  time-bandwidth  product  much  greater  than  one26 
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1 5  !(co)  |  =  const  =  , 

A 


and 


0(co)  =  const  =  8  . 


Substituting  (A.  18)  in  (A.  17)  yields 


r(r)  =  p  expO'0)  J  H (<W i)  exp(y' coz =  P  exp(j 0)  h (t  \t i)  . 


Substituting  (A.  19)  into  (A.  15)  gives 


where 


t/i(y)= -r  J  Aa51(05*(0^+  J  pexpO'0)/i(T;/i)5*(x)dT  (A.20) 


2  2 


-r  _r 

2  2 

=  yO AT  +  K(y)  , 


(A. 18) 


(A.  19) 


K(tr)=  Pexp0'8)  J  h(x-,tl)S\(z)dz  . 

1  L 

~  2 

—  T 

As  A  (t;f )  is  a  complex  Gaussian  process  and  S*(f)  is  only  a  phase  shift,  x(y)  is  also 
a  Gaussian  random  variable  with  zero  mean  and  variance 

I  T 
^  2  2 

Var(K)=±E  [k(t)k*(5)]  =  -|-  j  j  S(t)S*©  £[  h(v,t)  h'  (4;/)  ]  dx  d\  ,  (A.21) 

Z  Z  T  T 

2  2 

where  £[  •  ]  denotes  the  expectactation  of  [  •  ].  Assuming 

£  h(x\t)  /j*(^;f)  j  =  ct*25(t-£) 

we  have  Var(K)  =  al  T  . 


T  T 

It  is  obvious  that  if  k(— )  is  a  Gaussian  random  variable,  then  U  i(— )  is  also  a 


Gaussian  random  variable  but  with  mean  =yaAr  and  o2^  oj;T.  Similarily, 


we 


find  that 
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uA=\aA  j  S.iOSliOdt  +  j  P  expO  <3)  h(x\t)S*0(t)  dx  (A.22) 


=  yaA  J  cos(fir2)  dt  +  (3  exp  (7  0  )  J  cos(pT2)  h(x\t)dx  . 


j  1 

f/0(— )  also  is  a  Gaussian  random  variable  with  mean  =  -yot A  J  cos(|i/2)  dt  and: 


L 

R2  f 

a2  =  -^-a2  j  cos2(|i/2)  dt  .  (A. 23) 

2  r 

2 

T  T 

For  fixed  values  of  U0 (—)  and  t/i(— ),  the  probability  of  correct  decision  Pc  is  given 
by 


where 
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P  j  T  — 
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'Fh>'ffo 


‘A  exp  [  -V4  (Yo  +  Yi)l  / ot  'ftiy'o  1  , 


Yi  = 


p?<f> 
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(A. 24) 


Yo  s 


2  CT fj 


Then  the  marginal  probability  of  correct  decision  then  is: 


—  f  [  P  T  T 


P  T  dU i(4-)  P  . 

(tf.(y)  2 


L)dU  o(^) 


(A. 25) 
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A.4  Chirp  differential  BPSK  as  an  efficient  diversity  technique 

A  new  modulation  technique  now  will  be  presented.  This  technique  com¬ 
bines  the  chirp  waveform  and  the  phase  modulation  approach.  As  mentioned  previ¬ 
ously,  coherent  detection  in  the  presence  of  fading  is  useless.  This  especially  is  true 
when  the  data  is  transmitted  in  short  bursts  interspersed  with  long  breaks.  Even  if  the 
information  could  be  sent  coherently,  phase  synchronization  would  have  to  be  per¬ 
formed  on  each  burst.  Since  bursts  are  relatively  short,  it  is  difficult  to  achieve  phase 
synchronization. 

As  the  fading  is  very  slowly  time  varying  and  the  relative  drift  between  the 
transmitter  and  the  receiver  is  assumed  to  be  very  small,  the  phase  difference 
between  two  consecutive  bits  also  is  small.  Nevertheless,  it  must  be  taken  into  con¬ 
sideration. 

In  this  section,  we  will  prove  that  the  new  chirp  differential  BPSK  (CDPSK) 
is  identical  to  the  conventional  DPSK  in  the  presence  of  AWGN.  Then  we  will 
evaluate  the  performace  of  the  receiver  under  the  condition  of  fading.  Two  kinds  of 
fading  will  be  examined.  The  first  is  very  slowly  time-varying,  frequency  selective, 
Rayleigh  fading  and  the  second  is  non-random  fading. 

Let  us  define  our  signal  set  as 

5,(/)  =  sin(o)o  /  +  lAiu2  +  <fo(0  +  8]  z  =0,1  ,  (A.26) 

where  <J>,  =  y  when  a  mark  is  transmitted  and  <)>,  =  -y  when  a  space  is  transmitted. 

z 

e  is  random  phase  uniformly  distributed  between  0  and  2k  .  Figure  A.2  shows  the 
block  diagram  of  the  receiver.  The  received  signal  plus  noise  is  correlated  bit  by  bit 
with  a  one-bit  delayed  version  of  the  signal  plus  noise.  The  output  of  the  correlator  is 
compared  with  a  threshold  set  at  zero.  A  decision  is  made  in  favor  of  a  1  or  a  0 


depending  on  whether  the  correlator  output  is  positive  or  negative. 


Figure  A.2  Binary  chirp  DPSK  demodulator. 


In  our  analysis,  we  assume  that  a  mark  is  transmitted  and,  without  any  loss 
of  gennerality,  we  assume  that  0  =  0.  Hence,  the  received  signal  is 

ei(t)  =  ASi(t)  +  n(t)  ,  (A. 27) 

where 

5, (f )  =  sinfcoo  f  +  'A|i/2+y  ]  (A. 28) 

=  cos[co0r  +  ‘Apr2  ] 

and  n(i)  is  additive  narrow  band  Gaussian  noise.  In  complex  form  this  is  given  by 

Si(t)  =  A  Rej  .S(f)esp(/a>oO  j  .  (A. 29) 

Expressing  ex{t)  and  e2(t )  in  their  quadrature  components  yields 

e,(t)=  [A  cos(V*p/2)  +  nc//(r)]  cos(o)0/)-  A  sin(*/ip/2)  +  ns!i(t)  sin(WoO  ;  /=1,2(A.30) 
Following  Stein’s  proof,  81  we  find 

P,  =  'Aexp(-yb)  .  (A.31) 

We  see  that  in  the  absence  of  fading,  the  chirp  DBPSK  system  has  the  same  perfor¬ 
mance  as  the  conventional  DBPSK  system. 

For  a  fading  channel  we  substitute  r(t)  in  (A.31)  instead  of  S(t)  where  r(n  is 
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defined  in  (A.  16).  Again  following  Stein’s  proof81 

e, (0  =  t  n (f )  +  Xi ( t )]  cos(cdo  sin(a>o  t  )  i  =1 ,2  ,  (A.32) 

where  xt  (r )  and  y,  (> )  are  the  quadrature  components  of  the  noise.  The  fading  is 
assumed  to  be  a  slowly  time-varying,  frequency  selective  process,  so  at  each  fre¬ 
quency  the  fading  is  assumed  to  be  constant  for  the  relevant  time  period  (i.e. 

As  the  frequency  of  the  chirp  signal  is  swept,  the  amount  of  fading  varies 
accordingly.  This  actually  is  the  reason  why  the  amplitude  of  the  received  signal 
changes  during  a  oit  period.  Since  the  channel  is  very  slowly  time-varying,  the 
amplitudes  are  the  same  from  bit  to  bit,  i.e. 

r\(.t)~r2(t)  =  r(t)  .  (A.33) 

Eb 

Let  us  define  yb  as  y*  =  —  ,  where, 

No 

T 

2 

Eb  =  lh  J  r(t)r’{t)  dt  . 

zL 

2 

Assuming  no  intersymbol  interference 

Eb  =  '/2  j  (A. 34) 

+*•  -Mo 

=  j  j  H(f  -,t)S(f)H'(f;t)S'(f)df  =j  f  \H(f;t)\2  \S(f)\2 df  . 

For  sr»l. 


r  a2t 


Then 


\S(f)\2=l 


2 B 
0 


else 


(A.35) 
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B_ 

E„=^^-  \  \H(f;t)\2df  (A. 36) 

^  “  B 

2 

and 

A  A 

2  2 

=  J  I^C/JOI2^/  =  J  |//(/;0l2#  •  (A. 37) 

W0  2  BW0  Jg  B  j>g 

2  2 

Since  W(/;r)  is  a  zero  mean  Gaussian  random  process,  | H(f  ;r)l  is  a  Rayleigh  pro¬ 
cess  with  standard  deviation  2 aft . 

For  H(f\t),  with  pH  «  B  (  pH  defined  in  section  1.1  ),  y  can  be  approxi¬ 
mated  by 

7=^-2  I H (r^r+n p# ) | 2  p«  =  i  Y«  ,  (A.38) 

B  n= 0  2  n= 0 

where  K  is  defined  as  the  number  of  degrees  of  freedom  of  the  channel  and  is  equal 

D 

to  to  the  largest  integer  which  is  smaller  than  —  and: 

P  H 

Y„  =  f-pW  |//(-y+n*)l2=y  Pw|//„|2  . 

Sampling  \H(f  ;/)  I  at  rate  greater  than  or  equal  to  pH  guarantees  that  the  samples  will 
be  uncorrelated. 

Assuming  that  \Hn\  is  a  Rayleigh  distributed  random  variable  with  variance 
equal  to  2ct«  and  making  the  the  appropriate  transformation,  we  find  that 

p(7n)  =  —7—  exp 

JbpH 

For  a  given  sample  of  r(t ), 

P,^-  xh  exp  ( -7 )  (A. 40) 

and 


B  7, 


YfcPtf 


U{  7n) 


(A. 39) 


/*.=  J  P"rP(1)d  7  • 


(A. 41) 
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K- 1 

The  characteristic  function,  <t»(iw -y) ,  of  7  when  y  is  defined  as  7=  £  yn  ,  is 


K_x  r  pin  '  -k 

<&(iuy)  =  E  [exp(i'u7)]  =  {"[£  [exp(iit7*)]  =  - 7  iu  + - 7  .  (A.42) 

n=0  [  YfcP«<Jr  J  L  %PwCtr2_ 

The  probability  density  function  of  7  is  the  inverse  Fourier  transform  of  the  charac¬ 
teristic  function  4>(iuy),  and  is  given  by 


Ptf)  = 


YiPtfOr 


-J— exp( - Z-j)U  (7) 

(*-1)!  7iP«crr2 


(A.43) 


A  B  • 

Substituting  yb  =  yba},  and  K  =  —  in  (A.43)  yields 

Pw 


Pfi)=  7 ^7T  exp  ( -  41)  U(y)  .  (A. 44) 

.  y*  J  yb 

Figure  A.3  gives  the  calculated  density  function  of  7  for  various  values  of  K  and  for 
yb  =  10  dB .  Figure  A.4  compares  the  simulation  results  and  the  calculated  results  for 
K  =4  and  76  =  10 dB  . 


Probability  Density 


0.20 


0  10  20  30  40  50  60 


ElbiO/No 

Figure  A.3  The  calculated  probability  density  function  of  y  for  various  values  of  AT 

and  y4  =  10  dB . 
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Figure  A.4  The  probability  density  function  of  y  for  K  =  4  and  yb  =  10  dB  .  Comparision 
between  calculated  and  simulated  results. 


Substituting  (A.44)  in  (A.41),  changing  variables,  and  integrating  yields 


As  K  -» °°  (uncorrelated  process): 


(A.45) 


P.  =  * 


2  exP  Hfc  )  • 


(A. 46) 


Figure  A.5  gives  the  bit  error  probability  of  a  chirp  DPSK  as  a  function  of  y  for  vari¬ 
ous  values  of  K  . 
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Bit  Error  Probability  og  function  or  the  SNR 
K=2,  4 . 18,20 


E(bi1)/Nb  (dB) 

Figure  A.5  The  bit  error  probability  of  a  chirp  DPSK  as  a  function  of  y  for  various  values 
ofK  . 


For  a  given  signal  to  noise  ratio,  the  probability  of  a  binary  bit  error  for  4 
phase  DPSK  with  Gray  coding  is  given  by75 

P«/r  =  c[«>&]  exp  —  (a^b2)  ,  (A.47) 

where 


a  =  ypL  (V2+V2-V2^V2) 


and 


b  =  yjl  (V2+V2+V2^2) 
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After  some  algebraic  manipulation,  (A.47)  is  identical  to 

<A-48> 

The  marginal  bit  error  probability  P,  is  given  by  (A.41)  and  Py  is  given  by  (A.44). 
Substituting  (A.44)  and  (A.48)  into  (A.41)  and  integtating  with  respect  to  y  yields 


7 

'  4k  J- 


9=0 


2+-ir--  'ftjb  COS(9) 
7 b 


(2-V2cos(0))-!  d  9 


(A.49) 


and  as  K  -»°o, 


rr  LK 

=  f  (2-^cos(0)r1exp[-Yi(2-V2cos(0))]40  .  (A.50) 

9=0 

Figure  A. 6  gives  the  calculated  bit  error  probability  as  a  function  of  the  signal  to 
noise  ratio  for  various  values  of  K. 


Figure  A.6  The  bit  error  probability  of  a  4-Phase  chirp  DPSK  as  a  function  of  y 
for  various  values  of  K . 
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For  non-random  fading,  the  received  i-th  signal  is 

5,(r)  =  A  [1  sin(  0)0f  + + <!>,■(/)) +  asin  ( CO0F+ ^  tn2  +  <t>i(t)  + 9 )  ]  ,  (A. 51) 

where  T=t  +td  ,  and  td  is  the  time  delay  between  the  direct  and  indirect  received  sig¬ 
nals  a  is  the  attenuation  of  the  indirect  signal  and  9  is  the  phase  difference 

between  the  direct  and  indirect  signal  caused  by  reflection  and  the  difference  in  path 

length.  Assuming  td«T  ,  and  no  intersymbol  interference,  the  received  (i-l)-th  sig¬ 
nal  is  given  by 

)  =  A  [  sin  ( coo  t  +  Viji/2  +  ))  +  a  sin  ( co0  T  +  XA  pi2  +  )  +  9  )  ]  (A.52) 

and 

Si(t)  =  A  [  sin  (  coq  t  +  lA\xi2  +  4>i (f ))  +  a  sin  (  g>o  F+  A  (if2  +  (t )  +  9 )  ]  .  (A. 53) 

For  simplicity,  let  us  define  q(/)  =  lA\xt2  and  £(r)  =  a»o td  +  p  t  td  +  Vi\i  tj  +  9 . 

Without  any  loss  of  generality,  we  can  assume  that  i-1  =  1,  i-22,  and  <}>,  =  <t>2=y  • 
Rewriting  (A.52)  and  (A.53)  in  their  quadrature  form  yields: 

«i(0=  A  [cos(^(t))  +  acos(C(0]  +  «ei(0]  cos(coo(0)-  (A. 54) 

A  (sin(£(r ))  +  a  sin  (  CO ))]  +  ns ,(/ )  j  sin  ( coo(r )) 
and 

e 2(0  =  A  [  cos(£(f ))  +  a  cos  (  £(r )]  +  ne2(t )  j  cos  ( coo(r ))  - 
A  [sin(£(t))  +  a  sin  (  £(0)1  +  nj2(0  j  sin  ( <Bo(0) 

It  is  convenient  to  rewrite  (A. 54)  in  terms  of  a  another  nonzero  mean  Gaussian  pro¬ 
cess.  Then, 

e;  (0  =  (t ) cos ( %{t )) - y, (0  sin ( £(t ))  =  Re  z,exp(;£(f)  i=l,2  .  (A.55) 

Again  following  Stein, 81  we  find 


P,  =  Vi  exp  ( -y )  , 


(A. 56) 
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where  y=  „  and 
2 N0  T 


q  =  A 


r  os  r 

J 

-0.5T 


[  1  +  a2  +  2a  cos  (  C(f)] 


dt 


|  ,  1  ,  (if dT  1 

-  A  T  1  +  a2  +  '’a  co'CtOo/j  -1-  -p/i  +9)  sir.c  (—  )  j  . 


From  (A.57)  we  can  show  that 


(A. 57) 


1-lb 


,  l  ,  \ltdT 

1  +  or  +  2a  cos(a>0td  +  —lit/  +0)  sinc(— — ) 


For  a  very  short  time  delay  for  which  \itdT  is  very  small,  we  can  write 

1  2 


Y=Y b 


1  +  a  +  2a  cos(coofi  +  0) 


(A. 58) 


In  this  case,  y  is  pericdic  with  respect  to  td  with  a  period  equal  to  Tc  .  For  appropri¬ 
ate  values  of  td ,  coq,  and  0,  y  is  bounded  by  0  and  2yb,  and 


^<P.  <  y  exp  (-27*)  . 

In  most  underwater  acoustic  channels,  the  dominant  indirect  signal  is  the 
component  that  is  reflected  from  the  sea  surface.  In  this  case  0  =  k  .  Figure  A. 7  shows 
the  effect  of  the  indirect  path  on  the  bit  error  probability  for  various  values  of  time 
delay,  td,  between  the  direct  and  the  indirect  path  when  a  =  0.5  ,  and  0  =  n. 


Figure  A.7  Indirect  path  effect  on  the  BER  of  chirp  BDPSK  modulation  for  a  =  0.5 
and  0  =  n. 

fitjT 

Usually,  td  is  such  that  sinc(-^—)  =  0.  In  this  case: 

Y  =  Yi,(l  +  a2)2  (A.59) 

From  (A.59),  one  can  see  that  when  a  chirp  waveform  is  used  as  the  carrier  signal,  a 
short  delayed  deterministic  multipath  does  not  cause  any  dammage  to  the  system. 
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A.6.1  Motion  effect  on  the  performance  of  chirp  DPSK 


In  most  communication  links,  a  relative  motion  between  the  receiver  and  the 
transmitter  exists.  In  a  communication  link  between  several  sonobuoys  which  are  not 
moored,  currents  in  the  ocean  cause  a  relative  drift  among  the  sonobuoys.  Although 
the  time  drift  between  two  successive  bits  caused  by  the  distance  difference  is  very 
small,  it  is  a  mistake  to  ignore  this  effect  on  the  system  performance.  The  drift  is 
assumed  to  be  much  much  smaller  than  the  bit  period. 

Let  us  assume  that  the  i-th  received  bit  is  S,-(0  and  the  (i-l)-th  received  bit  is 
S (,_!)(/),  and  the  time  delay  caused  by  the  drift  is  much  smaller  than  the  bit  period. 
Then, 

5,(f)  =Ai(t)  sin  [  city  +  Vi  p/2  +  (0  +  0  ]  (A. 60) 

and 

S,_i(f)  =  /WO  sin  [  <%?  +  xh  |i/2  +  +  9  ]  . 

For  £  as  defined  in  (A.52)  and  5  =  coo td  +  \md  +lA  \itd  -r  e 
we  find  that 

S,(0  =>*,•(/)  sin  5(0  ,  (A. 61) 

and 


Si-iti )  =  A,_i(f )  sin  [5(f )  +  C(0 1  ■ 

From  Pawula  at  el70.  [(9)  -  (11)  ],  we  find  that  the  bit  error  probability  Pt 


is 


Pt  =p  (<!>,  <0<02)=^ 


f  (Qi)-F  (*,  )+i  ;  <  C(- j) < ^2 

F(<D2)-F  (<t>i)  ;<&,>C(f)or  <t2<C(y) 


where  <t>,  =  —  and 


A  3re  rc 

$7  =  —  = 


(A.62) 


and 
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K  T 

sin  [  £(— )  -  O  ]  f  exP  -YH-cos(C(j)-0))cos(f  ) 

F(<t>)= - JL - j  — L - - Ldt  {A63) 

—  1  -  cos  (  C(y )  -  <t>))  cos  ( o 


For  <j>;  =  0,_,,  a  correct  decision  is  made  if  <  d>  < 


Define  P/y  as 


P/y=F(-j)/y-F(f)/y  = 


(A. 64) 


cos(£(-J))  f  exP  _7  [1  _  sin  cos  (r )]  exp  -y[l  +  sin  («■£))  cos  (/)] 

2  (•  L  2 

-  - == - + - L _ _ 

4  TC  J  T  T 

-f  l-sin(C(j))cos(/)  1  +  sin  ( C(  j))  cos  (0 


The  marginal  probability  P  is 

^  =  J  P<yP(y)dy  (A. 65) 

where  /»(y)  is  given  by  (A.45).  After  some  simple  algebraic  manipulation,  one  finds: 


cos  (C(-jr))  2  ,  r  -  n 
f  i,  1  ,  .  7 bai 


!C  1=1  “I  * 


(A.  66) 


where 


ai  =  1  -  sin(C( y))cos(f) 


a2  =  1  +  sin(^(y))cos(f) 


COS(C(j))  2  2 

P  =  — 4^ —  /  I  ~ exp 

a  it  i=l  ut 


(A. 67) 
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From  (A.61),  we  find  that  the  bit  error  probability  Pe  is  given  by 


P.= 


f  >«{)«• 


(A.68) 


When  there  is  no  delay,  £(y)  =  0.  Substituting  this  into  (A.66)  and  (A.67) 

and  integrating  gives  us  (A.32).  Figure  A.8  shows  the  bit  error  probability  as  a  func¬ 
tion  of  the  signal  to  noise  ratio  (yb)  for  various  values  of  the  delay  7. 

drift  velocityi  0,  0.3,  0.6, . 1.5  m/sec 

Bit  Errcr  Probability  os  o  function  of  the  SNR  Jar  K=10 


Figure  A.8  Motion  effect  on  the  BER  of  chirp  BDPSK  modulation. 


Appendix  B  -  Performance  of  several  common  error 

correction  codes 

The  folowing  plots  taken  from  Odenwalder66  show  the  pperformance  of 
several  common  error  correction  codes. 


Figure  B.l  Bit  error  probability  versus  channel  error  probability  for  block 
length  n  -2m  - 1  Hamming  codes  with  m=3,4,  and  5 
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Channel  Error  Probability  ■  p 


Figure  B.2  Block,  and  bit  error  probabilities  versus  channel  error  rate  with 
extended  Golay  coding. 


Figure  B.3  Bit  error  probability  versus  channel  error  rate  performance 
of  several  block  length  127,  BCH  codes  C 
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Figure  B.4  Bit  error  probabilit 
of  several  convolut 


y  versus  channel  error  rate  performance 
ional  coding  systems. 


Appendix  C  -  Proof  of  Equation  3.16 


The  mean-square  error,  e\  between  the  transmitted  and  the  received  data  is 

e2  =  E^lV(n)  -  V(/.)J2j  =  £  j  [  £  (ctj  -tu)  2'  J2|  .  (C.1) 

Changing  order  of  summation  and  the  expectation  operation,  we  get 

e2  =  I  Z  E  f  (Cj)  -C'j )  (cnJ  -  emJ ) }  2f+"  .  (C.2) 

i  -O  ft  -0  L  J 

For  i  *  n 

E  [  (C,y  -*y)  (C„y  -  ^.y)]  =  £  [  (C,y  -*y)]  £  [  (C,„  -<■„)]  .  (C.3) 

But 

E  (Cjj  -Cjj)  j  =  0  . 

Since 

e[  (Cy-^y)]  =  j  [/»,  (1-0)  +  (1-/7,)  (1-1)  +/7,  (0-1) +  (!-/>,)  (0-0)]  =0, 

where  p ,  is  the  probability  that  the  i-th  bit  was  incorrectly  detected.  In  the  same  way, 
we  get  for  i  =  n , 

E  [  (c,y  -*y)1]  =Pi 

and  hence 

<‘■3-1(1?  •  (C.4) 

L  I*) 

Q.E.D 
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1 

1 

1 

Appendix  D  - 

Optimal  Power  Assignment 

1 

A  set  of  pi  is  to  be  found  such  that  the  function  e2  given  by 

1 

1 

1  4-1 

L  i-O 

1  b~l 

=  Z/(7,  )2'  . 

z  i-0 

(D.l) 

will  be  minimized  under  the  constraint  of  a  given  value  of  word  energy. 

Let  pi  be 

1 

the  probability  of  incorrect  detection  of  the  i-th  bit  in  the  word  (W).  It  is 

a  function 

i 

of  the  appropriate  signal  to  noise  ratio  y(.  In  other  words,  we  want  to  minimize 

1 

h  =  £/<*>  2* 

i*) 

(D.2) 

1 

1 

under  the  constraint 

4-1 

0=  E  y>  -y=o- 

f  -0 

(D.3) 

1 

Ew 

where  y  =  -f-  and  Ew  is  the  word  energy.  Using  Lagrange  multipliers,  we  define 

*M> 

< 

II 

J3 

1 

>* 

CD 

(D.4) 

1 

1 

1 

Then  the  set  of  y,  that  satisfies 

< 

0  =  0 

(D.5) 

■ 

1 

is  the  set  that  minimizes  (D.2).  In  our  case  , 

■ 

1 

$..?■  9r™  +x . 

(D.6) 

■ 

1 

1 

For  binary  DPSK  where 

/(yf)  =  y  exp(-y,)  . 

(D.7) 
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the  optimal 


Substiuting 


where: 


7;  are 


7 i  =  -£-  +  (2/  +  1  -b)  In  2, 
b 


7 >b(b  -  l)ln2 


(D.8)  in  (D.l)  we  get 


2  b  22-1  ,  „ 

rop,=  „h  ln2’ 


[2®  -  l]2 


7 >b(b-  l)ln2 


(D.8) 


(D.9) 
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